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1 Part 1: Singular points of solutions to analytic differential 
equations 



(i.i.i) 



1.1 Differential equations in the complex domain: basic results 

Recall that an analytic (also holomorphic) function w = u + iv of the complex variable 
z = x + iy is a differentiable map 

m 2 -> m 2 

(x,y) i ^ (u,w) 

satisfying the Cauchy - Riemann equations 

du dv du dv 
dx dy 1 dy dx 

Introducing the complex combinations 

z = x + iy, z = x — iy 
one can recast the Cauchy - Riemann equations into the form 

dw _ 1 / dw . dw\ 
dz 2 \dx dy J 

Theorem 0. Any function w(z) analytic on a neighborhood of a point z = zq admits an 
expansion into a power series 

w(z) = w + wi(z - Zq) + w 2 (z - z ) 2 + . . . (1.1.2) 

convergent for \z — zq\ < e for some sufficiently small e > 0. The coefficients of (1.1.2) can 
be computed via the derivatives of w(z) as follows: 

1 cftw(z) 

w = w{z ), w k = — — k | z ^ Zo , A; = 1,2,.... (1.1.3) 



Remark 1.1.1 The function w{z) is sais to be analytic at the point z = oo if the function 
w(l/u) is analytic at the point u = 0. Such a function can be expanded in a series of the 
form 

w = W0 + ^ + ^ + ... (1.1.4) 

z Z A 

converging for \z\ > R for some positive R. 

A function f(w, z) of two complex variables is called analytic if two Cauchy - Riemann 
equations in z and w hold true: 

df df 

-4 = 0, = 0. 

dz ' dw 
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Functions of two variables analytic on a polydisk 

\z\ < r, \w\ < p 

for some positive r, p admit expansion in a convergent double series 

f(w,z)=J2 a kiw k z l . (1.1.5) 

k,l>0 

Let us begin with considering a first order ODE of the form 

w' = f(w,z) (1.1.6) 

with a function f(w, z) analytic in a neighborhood of the point (wq, zq). We are looking for 
a solution w = w(z) satisfying the initial condition 

w(z ) = Wq. (1-1-7) 

Without loss of generality we assume that zq = wq = 0. 
Let us look for a solution in the form of power series 

w(z) = WQ + W\Z + W2Z 2 + . . . . (1.1.8) 

Substituting into the equation (1.1.6) one obtains a recursive procedure uniquely determining 
the coefficients: 

w = 

Wl = w '(0) = /(0,0) = aoo 

W2 = \ ( ) = \ [f z ( w > ^ + f w ( w ' z ^ w '^ (o,o) = \ ( fl01 + a io°oo) 

w 3 = \\-r [fz( w ' z ) + fw{w,z)w']) 

6\dz L J 7(0,0) 



/z^ + + w' 2 f ww + w"f, 



(0,0) 



= g [^02 + 2ana o + 2a 20 aoo + «io(«oi + aioa o)] • (1.1.9) 

In these computation Ojj are the coefficients of the Taylor expansion (1.1.5) of the function 
/■ 

In general 

w k = P k ( aij ) (1.1.10) 

where Pk is some polynomial with positive coefficients. It remains to prove convergence of 
the series. 

Exercise 1.1.2 Prove that the series 

w = z + z 2 + 2! z 3 + 3! z 4 + • • • + (n - 1)! z n + . . . (1.1.11) 
satisfies the differential equation 

z 2 w' = w — z. 
Prove that the series diverges for any z ^ 0. 



The following classical result, due to Cauchy, establishes convergence of the series solution 
(1.1.8), (1.1.9) for differential equations with analytic right hand side. 

Theorem 1.1.3 // the function f(w,z) is analytic for \z\ < r, \w\ < p for some r > 0, 
p > then the series (1.1.8), (1.1.9) converges to an analytic solution to (1.1.6) satisfying 
the initial condition (1.1.7) on the domain 

\z\ <n(l- e~^n) (1.1.12) 

for arbitrary r±, p\, M such that 

<n < r, < pi < p 

assuming that 

\f(w,z)\ < M for \w\ < pi, \z\ < n. 

Let us begin the proof with a definition. Given two power series 

f(w, z) = J2 a kl w k z l , F(w, z) = Y J A M w k z l (1.1.13) 

such that (i) the series F(w,z) converges for \z\ < n, \w\ < pi, and (ii) all the coefficients 
A^i are real positive numbers, we say that F(w, z) is a majorant for f(w, z) if 

\&kl\ < Aj~i for all k, I > 0. 

We will use the following notation 

f(w,z) ^ F(w,z) (1.1.14) 

to say that F(w, z) is a majorant of f(w, z). Observe that the power series f(w, z) converges 
on the same domain \z\ < n, \w\ < pi where the series F(w,z) does converge. 

Lemma 1.1.4 Given an arbitrary polynomial P(aoo, aoi, • • • , a^i) with positive coefficients, 
then for any pair f(w, z) ■< F(w, z) the following inequality holds true 

|P(a 00 , a i, • • • , a k i)\ < P{A 00 , A 0h . . . , A kl ). (1.1.15) 

The proof is obvious. 

Let us now recall the Cauchy inequalities for the Taylor coefficients of an analytic function. 

Lemma 1.1.5 Given a function f(w, z) analytic for \w\ < p, \z\ < r, then for any < r\ < 
r, < pi < p the following equality holds true 

-| r2ir p2tt 

Y,Wki\ 2 pl k rf = -—^ / \f{p 1 e i tr l e i6 )\ 2 d ( j>de. (1.1.16) 
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Proof: Indeed, 



/(pie* ne*) =J2^W k Ae- 



ik<f>—iW 



Multiplying one obtains 
\u>kl\ Pi r 21 

k,l 

where we use the orthogonality 

r-2ir flit 



JO JO 



□ 



Corollary 1.1.6 Denote 

M= sup |/(w,2)|. 

M<Pl> |z|<ri 

T/ien i/ie following inequalities hold true for the Taylor coefficients of the function f(w, z) 

M<^j Vfe, />0. (1.1.17) 



Proof: Indeed, 

]T|a H | 2 pfrf<M 2 . 

Hence 

□ 



Corollary 1.1.7 

M 

f(w, z) < V7 v . (1.1.18) 



Proof: follows from 

^pM 1 "" "(i-^) (i-f) 



2 y -TTTW 2 - 



□ 
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Proof of Cauchy theorem: consider an auxiliary differential equation 

w. M 



1 



1 



111 

pi 



(1.1.19) 



Let us consider the power series solution to this equation 

W = B mZ rn , W(0) = 0. 

m>l 

The coefficients B m can be determined by the procedure (1.1.9) in the form 

B m = P m (boo, boi, ■ ■ ■ ) 
by the same polynomials P m as in (1.1.10). Here 

7 M 

hi 



p\r\ 



are the Taylor coefficients of the rhs of (1.1.19). According to Lemma 1.1.4 

w{z) < W{Z). 

It remains to prove convergence of the power series solution to (1.1.19). This equation can 
be solved explicitly by separation of variables: 



Pi J 



1 - 



This gives 



1 - 



W 
Pi 



2Mr ± 
pi 



log 1 - 



Z 

n 



+ K. 



The initial condition W(0) = determines the integration constant 

K = 1. 

Finally 

W = pi 



This function has singularities at Z = n and 

Z 



2Mn 



pi 



log 1 - 



n 



(1.1.20) 



log 1 - 



Pi 
2Mri 



i.e. at 



Z = n(l-e 2 ^^i) 



(observe that < Zo < n). Hence the Taylor series for the function (1.1.20) converges for 

\z\ <n(i- e _5 ^^r^) . 

Due to w(z) < W(Z) the series solution w(z) converges on the same domain. □ 
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Remark 1.1.8 Uniqueness of the analytic solution easily follows since the procedure (1.1.9) 
determines all the coefficients uniquely. 



Remark 1.1.9 Let the rhs of the differential equation 

w = f{w,z) 

be defined on the (punctured) polydisk 

\w\ < p, \z\ > R 

for some positive p, R. We say that the point z = oo is a point of analyticity for the 
differential equation if, after the change of independent variable 

1 2 du> t> / -i / \ 

z = -, -u — = f (w, 1 u) 
u du 

the point u = is a point of analyticity, i.e., the function 

f(w,u) := -u~ 2 f{w,vT l ) 

is analytic for 

1 

H < P, \ u \ < ^- 

In that case the Cauchy theorem ensures existence, uniqueness, and analyticity of solutions 
to the differential equation with the initial data 

w(oo) = Wq. 

The solutions are represented as power series in 1/z: 

w = w -\ 1 — =- + ... 

z z z 

convergent for \z\ > R\ for some positive R±. 

I will leave as an exercise to formulate and prove an analogue of the Cauchy theorem for 
systems of differential equations 

w[ = fi(w!, ...,w n ,z) 

(1-1-21) 

w' n = fn(m, ...,W n ,z) 

with analytic right hand sides. 

The Cauchy theorem can be improved in certain particular cases. As an important ex- 
ample consider a system of linear differential equations 

dw n 

= J2 a kl(z)wi + b k (z), k = l,...,n (1.1.22) 

l=i 

with coefficients aki(z), bk(z) analytic for \z — zq\ < r for some positive r. 

Exercise 1.1.10 Prove that the solution to the linear system (1.1.22) with the initial data 

wi{z Q ) = Wi,..., w n (z ) = w Q n 
with arbitrary (m®, . . . , w„) G C n exists, is unique, and analytic on the same disk \z — zq\ < r. 
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1.2 Analytic continuation of solutions. Movable and fixed singularities. 

The considerations of the previous section describe the local structure of solutions to differ- 
ential equations near the point of analyticity. Global considerations require working with 
multivalued solutions to differential equations that we will introduce in this section. 

We begin with recalling basic facts about analytic continuation. 

Definition 1.2.1 An element of analytic function is a pair (U, /) where U C C is a disc and 
f is an analytic function in U. 

Definition 1.2.2 An element (V,g) is called analytic continuation of (U, /) if 

1) unv / 

2) on the intersection U fl V the identity g(z) = f(z) holds. 

Definition 1.2.3 A multivalued analytic function w(z) on a connected domain f2 C C is a 
set of elements of analytic functions 

(Ui,w id ) ieI; jeJ (1.2.1) 

such that 

1) UieiUi = V. 

2) For any pair Ui 17 Ui 2 with a non-empty intersection and for any j± G J there exists 
22 G J such that (Ui 2 ,Wi 2 j 2 ) is analytic continuation of (Ui 1 ,Wi 1 j 1 ). 

3) On triple intersections ofU^, Ui 2 , Ui 3 the analytic continuations Wi 2 j 2 andwi 3 j 3 must 
coincide. 

For a multivalued function on Q C C it is defined an operation of analytic continuation 
of an element (U^ , u>i 1 j 1 ) from a point P G Ui 1 to another point Q along a curve C. Namely, 
choose a finite number of elements (Ui 1 , w^j^, . . . , (Ui n ,Wi nt j n ) such that 

n 

Cc(JU ik , 
k=l 

PeU h , QeU in 

and the element (Ui k+1 , Wi k+1 ,j k+1 ) is an analytic continuation of an element (Ui k ,Wi k j k ) 

Exercise 1.2.4 Prove that the value Wi n j n (Q) of the analytic continuation of an element 
(C^ii, of a multivalued function from P £ to Q G Ui n along C does not change with 

continuous deformations of the curve C . 
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In particular, any closed curve 



C : {z = z(t), < t < 1, z(0) = z(t)} 

induces a bijection 

HC--J->J (1-2-2) 

denned by the following condition: for any j € J the element {Ui n i w i n ,jx c (j)) ls ^ ne result 
of analytic continuation of the element (f7j , Wj j) from the point P = zq to the same point 
P = z\ along the closed curve C. The bijection (1.2.2) depends only on the homotopy class 

[C]€7Tl(fi,P) 

of the closed curve C. In this way one obtains a homomorphism 

/i : 7Ti(fi,P) -> Aut(J) 

(1.2.3) 

[C] 

(1.2.4) 



Definition 1.2.5 T/ie homomorphism (1.2.3) is called monodromy o/ i/ie multivalued ana- 
lytic function. 

In particular, if the homomorphism fi is trivial, 

/i(7ri(n,P)) = id G Attt(J) 

then the multivalued function is actually a collection of J functions analytic in 17. 
An equivalent version of Definition 1.2.3 is 

Definition 1.2.6 A multivalued analytic function on O C C is an analytic function on some 
covering Cl of the domain 17 

Remark 1.2.7 The covering 17 has a natural structure of complex variety (of complex di- 
mension one). 

Back to differential equations: let us consider a more general differential equation 

F(z,w,w') = (1.2.5) 

where 

F(z,w,w') = Y J a k,i( z ) wk (w'Y 

k,i 

is a polynomial in w, w' with coefficients aki(z) analytic in z on a domain D G C. 
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Theorem 1.2.8 If w = w(z) is a solution to the differential equation (1.2.5) analytic on the 
disk U C D and (V,w) with V C D is an analytic continuation of the element (U,w) then 
also the function w(z) satisfies (1.2.5). 



Proof: Indeed, 

F(z,w(z),w(z)) = for z G U l~l V. 

Due to the "principle of analytic continuation" this function vanishes identically for all z <E V. 

□ 

The domain of analyticity and behaviour of solutions to analytic differential equations 
are determined by their singularities. The study of positions of the singularities and their 
properties is one of the main tasks of the analytic theory of differential equations. 

Recall that a point zq G C is an isolated singularity of a function w(z) analytic in if 
there exists a disk U such that 

u \ z c n 

and zq fi. 

Let us first subdivide isolated singularities into two big classes. 

Definition 1.2.9 We say that zq is a critical singularity of a multivalued analytic function 
w(z) if the analytic continuation along a nontrivial closed loop inU\zQ changes the value of 
the function, and non-critical in the opposite case. 

For example, z = is a critical singularity for functions *Jz, logz, but it is a non-critical 
singularity for the functions 1/z, l/z k . 

We will now consider isolated singularities of analytic differential equations. 

Definition 1.2.10 Let w(z) be a solution to the differential equation 

w = f(w,z) 

with an isolated singularity at the point z = z$. We call this singularity fixed if its position 
does not depend on the choice of the solution w(z); in the opposite case the singularity is 
called movable. 

Example. Consider the differential equation 



2w z 

The solution can be found explicitly: 




where C / is an integration constant. It has fixed transcendent singularities at z = and 
z = oo. It also has movable critical singularity at z = C. 

Observe that the fixed singularities z = and z = oo can be determined looking at 
the singular points of the coefficients of the differential equation. Localization of movable 
singularities is often a hard task. The situation simplifies for linear differential equations. 



10 



Theorem 1.2.11 Solutions of linear differential equations with rational coefficients have no 
movable singularities. 



Proof: Let z±, . . . , zn £ C be all the poles of the coefficients of the system. Any solution 
defined in a neighborhood of a point zq G C\ {z±, . . . , zjy} can be analytically continued along 
any curve on the punctured sphere. In this way we obtain a multivalued analytic function 
defined on all C \ {z\, . . . , zn}- □ 



Corollary 1.2.12 All movable singularities of Riccati equation 

w' = a(z)w 2 + b(z)w + c(z) (1.2.6) 
with arbitrary rational coefficients are poles. 

Proof: Consider the following linear differential equation 

y" + p(z)y' + q(z)y = (1.2.7) 

with 

p(z) = -b(z)- a ^, q(z) = a(z)c(z). 

The substitution f 

w{z) = -^ rs V - (1.2.8) 
a{z) y 

into (1.2.7) yields (1.2.6). Since the solutions y{z) to (1.2.7) have no singularities away 
from the poles of rational functions p(z), q{z) (these are fixed singularities), the movable 
singularities of w(z) are poles (at the zeroes of y(z)). □ 

We want now to prove the following converse statement. 

Theorem 1.2.13 A differential equation of the form 

w = f(w,z) 

with rational rhs 

_ a (z)w n + a 1 (z)w n ~ 1 j h a n {z) _ P n (w, z) 

J[W > Z) - (3 (z)w™ + (3 1 (z)w^ + --- + p m (z) - Q rn (w,z) [ ■ ■ > 

has no movable critical singularities iff it is a Riccati equation. 

Let us begin the proof from the following lemma (useful by itself) describing the local 
structure of solutions to analytic differential equations near the poles of the rhs. 
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Lemma 1.2.14 Let the rhs f(w,z) of the differential equation w' = f(w,z) have a pole at 
the point (wq,zq), i.e., the function l/f(w,z) is analytic near (wq,zq) and 



lim 

w — > u>o 



f{w,z) 



0. 



Then the differential equation has a solution in the form of a Puiseaux series 

w = w + ai(z - z ) 1/k + a 2 (z - z ) 2/k + ... 
for some integer k > 1 convergent for \z — zq\ < r for some positive r. 

Proof: Rewrite the differential equation for the inverse function 



dz 



1 



dw f(w,z) 

Applying Cauchy theorem one obtains the solution with the initial data 

z(wq) = ZQ 

in the form 



Observe that 



Denote 



z = z + b\(w - w ) + b 2 {w - wq) 2 + 

1 



bi 



f(w,Zjy {wihZo) 

k := min{i | bi / 0}. 



0. 



Rewriting the equation 



in the form 



z - zo = b k (w - w ) k + b k+1 (w - w ) k+1 + . . . 
z - zo = b k (w - wo) k [l + ci(w - w ) + c 2 (w - w ) 2 + • . . 

h+i , A . b k+ - 



1 + —r (W — W ) + — [W - Wo) + 

bk b k 



where 

l + c\{w — wo) + c 2 {w - wq) 2 H 

we apply the analytic version of the implicit function theorem to the equation 
C = h(w - wq) [l + a(w- wq) + c 2 (w - wq) 2 + ■■•], h := b]/ k 
One obtains an analytic function w = w(() for sufficiently small \(\: 

w = wo + ai( + a 2 C +■■-, a\ = — ■ 

bi 



(1.2.10) 



-i i/fc 
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The substitution 

gives the needed solution w(z). 



(=(z-z ) 1 / k 



□ 



We see that poles of the rhs of differential equations correspond to algebraic critical 
singularities of solutions. 

Example. The solutions to differential equation 



(n + l)w n 

have the form 

i 

w(z) = (z — Zq) "+ 1 . 

For any n > it has a movable algebraic critical singularity. 

Let us return to the proof of Painleve theorem. If the degree m of the denominator is 
positive then, for any zq G C such that (3o{zo) 7^ there exists a point (wo, zq) such that 

<5m(^o,^o) = 0. 
Without loss of generality one can assume that 

P n (w ,z ) 

(the polynomials P n (w,z) and Q m (w,z) by assumption have no common factors) and 

d w Qm(w,z)\^ WOjZ) / 0. 

The last condition ensures local analytic dependence on zq of the root w = wq of the poly- 
nomial equation 

Q m (w,Z ) = 0. 

So, according to Lemma 1.2.14 there exists a solution to the differential equation with a 
critical singularity at z = zq. The possibility of small variations of zq proves that this 
singularity is movable. 

We proved that m = 0, i.e., the differential equation must have the form 
w' = a (z)w n + ai{z)w n ~ l H h a n (z). 

Put 

1 

w = —. 

w 

One obtains 

_ ao(z) + a>i(z)w H h a n (z)w n 

yjn-2 

If n — 2 > then this equation has movable critical singularities, i.e., the original equation 
has movable critical poles. Hence n = 2 and we obtain a Riccati equation. □ 
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1.3 Singularities for differential equations not resolved with respect to the 
derivative. Weierstrass elliptic function. 

We consider only a particular class of such equations having the form 

F(w,w') = (1.3.1) 

for a polynomial F(w,w'). Let us begin with 

Example. Consider one-dimensional motion of a point of mass 1 in the field of a poly- 
nomial potential V(x). Recall that the equations of the motion have the form 

Conservation of the total energy reduces this second order ODE to the first one 

^± 2 + V{x)=E (1.3.3) 

where E is the integration constant. 

Question: when solutions x = x(t) to the differential equation (1.3.3) have no movable 
critical singularities in the complex t-plane? 



Lemma 1.3.1 The solutions to (1.3.3) have no movable critical singularities only if degV(x) < 
4. 

Proof: The substitution 

1 dx 1 dx 

x 1 dt x 2 dt 

gives 

K§)' = ** [*-"<*"'>] = a < 1JM > 

If the degree of the polynomial V{x) is greater than 4 then the rhs of (1.3.4) has a positive 
power of x in the denominator. 

One arrives at the following question: given a differential equation of the form 

w = , ngZ 

w n 

with a polynomial 

P(w) = a + aiw + . . . , 
prove existence of movable critical points under the assumption n > 0. 
To do this we separate the variables 

w^dw 

— , = dz. 
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Let us consider the solution with the initial data w(zq) = with arbitrary zq. It can be found 
in the form 

(a J 2 w 2 l + Cl u> + c 2 w + . . . = z - z 

n + 2 1 J 

where the (convergent) series in the square brackets is defined by the following formula 

1 + Cl w + c 2 w 2 + ■■■= (a Q 1 P(w))~ 1/2 . 

Inverting one obtains the solution w = w(z) with a movable algebraic critical singularity (if 
n > 0) 

2 3 

w = bo(z — zq) n + 2 + b±(z — zq) n + 2 + . . . 
for some coefficients bo, b±, . . . . □ 

We will now prove that for n < 4 all movable singularities of solutions to the equation 

w ' 2 = ao w n + aiw n ~ l + ■ ■ ■ + a n (1.3.5) 

are poles. The statement is trivial for n = 0, so we start from 
Case 1: n = 1. General solution to 

w = z — a 

reads 

\ 2 

, Z- Zq 

w = a + 



Case 2: n = 2. To solve the differential equation 

w' 2 = (w — a)(w — b) 

with a ^ b let us consider the algebraic curve 

p 2 = («; — a)(w — b). 

The curve admits a rational parametrization 

a - bs 2 
1-s 2 

{a-b)s 
1-s 2 ' 



if = 



P 



We have 



w — a 



so 



i.e. 
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Integration of this equation gives 

1 - e z ~ z ° 



1 + e z - z o 

Therefore the solutions w = w(z) are represented in the form 

w = R( e z ~ Z0 ) 

for some rational function R. 

Case 3. n = 3. Doing if necessary an affine change of dependent variable 

w i ^ aw + b 

we reduce the differential equation (1.3.5) to the form 

/2 3 

w = 4w — g 2 w — 53. (1.3.6) 

Here g 2 , g% are some complex number. We will construct a particular solution w = p(z) 
to this differential equation called Weierstrass elliptic function. We prove that this solution 
is meromorphic on the entire complex plane. All other solutions are obtained by a shift of 
argument 

w = p(z - zo), z £ C. 
Let us fix a pair of complex numbers u\ , UJ2 such that 

Im— >0. (1.3.7) 

Denote 

A := {z = 2muJi + 2nuj 2 \ m, n G Z} (1.3.8) 

the lattice on the plane generated by the vectors 2uji and 2^2 • Elliptic functions associated 
with the lattice A by definition are meromorphic functions on the complex torus 

T W1)W2 := C/A. (1.3.9) 

Spelling this definition out we identify elliptic functions with doubly periodic meromorphic 
functions on the complex plane: 

f(z + 2u 1 ) = f(z), f{z + 2uj 2 ) =f(z). (1.3.10) 

Exercise 1.3.2 Prove that any elliptic function holomorphic on C must be a constant. 
The Weierstrass elliptic function is defined by the following infinite sum 



p(z;ui,u 2 ) := -2 + X] 



Z' 

u>eA\o 



(Z — Uj) 2 LO 2 



(1.3.11) 



Theorem 1.3.3 The series (1.3.11) converges absolutely and uniformly on any compact in 
C \ A. The sum p(z\u\,U2) is a doubly periodic analytic function in z G C \ A with double 
poles at the lattice points. 
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Proof: We first prove the following 



Lemma 1.3.4 The series 

y — 

converges. 

The notation 

ujeA 

will be used for summation over all nonzero elements of the lattice A. 

Proof: For any positive integer n let us evaluate the contribution of the lattice points 
belonging to the boundary of the parallelogram 

Il n := {z = 2uisi + 2w 2 s 2 | - n < si, s 2 < n}. 

Denote d the minimal distance from the origin to the boundary points of the parallelogram 
IIi. Then the distance from the origin of any lattice point on the boundary <9II„ is greater 
or equal to dn, that is 



\co\ 3 d 3 n 3 
We have 8n lattice points on dU n , hence 

' 1 8 
^ IcjI 3 — d s n 2 

The estimate for the sum of the series readily follows 

V'— - - V — - — 

2-f | w |3 -«j3L n 2- 3d 3' 
wGA n=l 



□ 



We are now in a position to prove convergence of the series (1.3.11) on any disk \z\ < r 
away from the lattice points. For a given r > there is only a finite number of lattice points 
inside the disk of radius 2r. Considering only the points uj G A such that \u\ > 2r the 
following inequalities hold true 



1 1 




2uj z — 


z 2 


(z — uj) 2 UJ 2 




uj 2 {oj — 


z) 2 



\z(2 - z/uj) 



to 



3 " M 3 



|o;| 3 |l — z/uj\ 2 
So the convergence readily follows from the Lemma. 



< r-^r for 



\z\ < r, z $l A. 



17 



For a z inside the disk \z\ < r, z A rewrite the series for the p-function as follows: 



(*;fi,w2) = i+ J2 



weA, M>2r 



ujSA, 0<M<2r 
1 1 " 



(z — uj) 2 UJ 2 



(z — uj) 2 UJ 2 



The second series is a function holomorphic for \z\ < r. The first (finite) sum, along with the 
1/z 2 term, is a rational function, hence it is meromorphic on the disk \z\ < r. This prove the 
statement of the Theorem about the poles. 

Let us consider now the derivative 

p'(z-ui,u 2 ) = -2 V - — —3. 

^— ' [z — UJ) 6 

It is clear that this meromorphic function is doubly periodic: the shift z 1— > z + 2ku\ + 2Zo>2 
can be absorbed by a change of the summation indices 

mh^m + k, tjHn+I. 

In order to prove double periodicity of the Weierstrass function let us first observe that 
the difference 

p(z + 2kui + 21uj 2 ) - p{z) 

is a constant for any k, I G Z. Indeed, the derivative of this function is identically equal to 0. 
To compute the value of the constant p(z + 2ui) — p(z) it suffices to set z = —uj\ and then 
use that p(z) is an even function, hence 

In a similar way one proves that p(z + 2^2) — p(z) = 0. □ 



Let us now prove that the Weierstrass function satisfies the differential equation (1.3.6) 
with 

^ = 60 E'^4> 53 = 140^'^. (1.3.12) 



Theorem 1.3.5 The Weierstrass function p(z; uj±, 0J2) satisfies the differential equation 



(1.3.13) 



Proof: Using the geometric series 



1 1 2z 3z 2 4z 3 5z 4 

l + — + ^r + ^ + —r + ... 



(1-5)' 



UJ UJ 2 UJ 3 UJ 4 
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we derive the Laurent expansion of the Weierstrass function at z = 0: 



P[ ' z 2 + 20 + 28 ' 



Hence 



Therefore 

This implies that the function 



/ // n\2 4 2g 2 4g 3 

(p(*)) =?-5^--r + --- 

S/ \ 1 3(72 3 



/ 2 -4p 3 = -^§- 53 + 0(^ 
z z 



p' -Ap 3 + g 2 p + g 3 = 0{z 2 ) 

is holomorphic at z = 0. Due to double periodicity it is holomorphic everywhere. Thus it 
must be equal to a constant. Since the value of this functions at z = is equal to zero, the 
constant is equal to 0. □ 

It turns out that, choosing in a suitable way the periods u\ j2 of the lattice one can 
obtain a solution to the Weierstrass differential equation (1.3.13) with arbitrary values of the 
parameter g 2 , g 3 satisfying 

g 3 -27g 2 ^0. (1-3-14) 

Actually the generators of the lattice are defined as the periods of a holomorphic differential 
over a suitably chosen basis of cycles a±, a 2 & B.\{C\ r E) on the elliptic curve 

C : p 2 = 4w 3 - g 2 w - g 3 , (1.3.15) 

/ dw / dw 

2ui = (b — , 2uj 2 = <p — . 

J ai P Ja 2 P 

The details of this construction go beyond the scope of this course. 

For g\ — 27g 2 = the polynomial 4w 3 — g 2 — g 3 has a multiple root. The curve (1.3.15) 
becomes rational: 

p 2 = 4(io - a) 2 {w - b) 
for some a, b. The rational parametrization of the curve reads 

w =b + s 2 

p =2s(b-a + s 2 ). 

The corresponding differential equation w' 2 = 4(u> — a) 2 (w — b) integrates in elementary 
functions. 

Case 4. n = 4. In order to solve the differential equation 

w' 2 = qqw 4 + Aa\w 3 + 6a 2 w 2 + Aa 3 w + 0,4 = P±{w) 
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let us again consider the algebraic curve 



The substitution 

yields 
with 



P 2 = Pa{w). 



w = s + a 



P 



(1.3.16) 



p 2 w 4 = b^w 4 + b%w 6 + b2W 1 + b\w 1 + P^a 



,7.-2 , L ~.-l 



& fc = -^P 4 (fc) (a), fc = l,-..,4. 



If a is a root of the polynomial P4(«;) then one obtains a cubic 



One more substitution 
with 



p = 64 + b^w + b2W + biw" 



p = Ay, w = A x + /i 



01 3oi 



(1.3.17) 



reduces equation (1.3.17) to the Weierstrass normal form 



y = 4a; 3 - 52a; - g 3 



(1.3.18) 



with 



92 = ^7,(^-36163) = a a4+3al-4aia 3 , g 3 = -t^(96i6 2 &3 — 26|— 276f & 4 ) 



ao a\ CL2 
ai a 2 a 3 
a 2 a 3 a 4 



If the quartic polynomial P&{w) had no multiple roots so does the cubic polynomial in the 
rhs of (1.3.18). Therefore the curve (1.3.18) admits uniformization by Weierstrass elliptic 
functions: 

w =Xp(s)+fi 

(1.3.19) 

P =Ap'(s). 



As 



dw dw 
az = — = — = ds 
p p 



we finally obtain the general solution in the form 



w = a + 



Xp(z - zq) + //' 

It is clearly a meromorphic function in z. 

More general result is given by the following theorem of Hermite (1873). 
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Theorem 1.3.6 If the polynomial differential equation 

F(w,w') = 

has no movable critical points then the genus of the Riemann surface 

F(w,p) = 

is less or equal to 1. The solutions can be of the following three types: (i) rational functions 
in z; (ii) rational expressions in exponential functions; (Hi) rational expressions in elliptic 
functions. 

Sketch of the proof On the curve F(w,p) = consider the abelian differential 

dw ( a \ 
to = — (= dz). 

P 

Since for any abelian differential 

#(zeroes) — #(poles) = 2g — 2 

the differential has a zero (wo,po) for g > 1. Choosing a complex coordinate r near such a 
point 

W = W + W\T + W2T 2 + . . . 
P = PO + P\T + P2T 2 + ■■■ 

one has 

lo = [a k r k + a k+ iT k+1 + . . . ) dr 
for some positive k G Z, a k / 0. So 

dz 



and 



(h a k T k + a k+l T k+1 + . 



fc+1 



or 

/ (fc + l)(z-zo) \ 5Ti 
We obtain a movable critical singularity: 



/(A; + l)(z - z )\ fc +! 
w = w + wi f ^ -) +0[(z- z )fe+ij 



□ 



As an application of the Hermite's theorem consider the problem of finding meromorphic 
functions possessing algebraic addition theorem. We say that the function f(z) satisfies an 
algebraic addition theorem if there exists an algebraic function P(u, v, w) such that for all 
x, y £ C the following identity holds true 

P(f(x),f(y)J(x + y)) = 0. (1.3.20) 
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Corollary 1.3.7 Meromorphic functions f(z) possessing algebraic addition theorem are ra- 
tional functions in z, or in e az , or elliptic functions. 

Proof: Differentiating (1.3.20) in x and y obtain 

Puf'(x) + P w f'(x + y) = 
Pvf'{y) + P w f'{x + y) = 0. 

Hence 

P*(f(x), f(y),f(x + y)) f'(x) - P v (f(x)J(y),f(x + y)) f'(y) = 0. (1.3.21) 

Eliminating f(x + y) from the equations (1.3.20), (1.3.21) we obtain an equation of the form 

<f>(f(x),f'(x),f(y),f'(y)) = 0. 

Fixing y one arrives at a differential equation of above form by assumption having mero- 
morphic solutions. Due to Hermite theorem the solutions must be rational, trigonometric or 
elliptic. □ 

Example (L.Fuchs). Consider the differential equation 

w' 3 - 3w' 2 - 9w 4 - 12u> 2 = 0. (1.3.22) 

The associated algebraic curve reads 

F(p, w)=p 3 - 3p 2 - 9w 4 - 12w 2 = 0. (1.3.23) 

This irreducible equation is cubic in p. So the associated Riemann surface has 3 sheets over 
the complex u>-plane. 

Step 1: look for the ramification points for the projection (p,w) i— > w of the Riemann 
surface (2.1.19). To this end we have to first solve the system 



F(p,w) =0 

F p (p,w) =0. 
The second equations gives 

3p 2 — 6p = 0, hence p = or p = 2. 



(1.3.24) 



One obtains 



2i 

p = p\ = =>■ w = wi = or w = u>2,3 = 

v3 



(1.3.25) 



p = P2 = 2 =>- W = W4£ = ±i 

Let us study the structure of the Riemann surface near the points (pi,Wj). 
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1) Near (p±,wi) = (0,0) one has, at the leading order 

p ~ ±2 i w. 

Higher order terms are determined uniquely for both branches: 

p±(w) = ±2iw- \w 2 + 0{w 3 ). (1.3.26) 

One obtains two analytic functions p = p±(w). Hence (pi,w\) is a double point of the 
algebraic curve (2.1.19) (i.e., not a branch point). It is obtained by identification of the two 
points (1.3.26) of the Riemann surface. 

2) Near (pi,^) or (pi,W3) one can use the coordinate p as the local parameter: 



w = ± 



2i i y/3 2 ^3 



+ 0(p l ) uca. (p.w) - (0,±^) • (1-3.27) 



So, the points (pi,^) an d {pi,w^) are both second order branch points on the Riemann 
surface. 

3) Near each of the points (^2,^4) or (^2,^5) both branches of the Riemann surface 
(2.1.19) are holomorphic: 



or 

w 



4) At infinity one has to do the substitution 

1 1 

w = -, p = - 
u q 

to arrive at 

9^ 3 - u 4 + 12uV + 3gu 4 = 0. (1.3.30) 
Near u = the curve can be described by a Puiseaux series 

q = 3- 2 / 3 eu 4 / 3 - + O (u 11 / 3 ) . (1.3.31) 

Thus at infinity the curve (2.1.19) has a third order branch point. 

We obtain on the Riemann surface 3 branch points of orders 2, 2 and 3. The total 
multiplicity of the branching locus by definition is equal to 

b = (2-1) + (2-1) + (3-1) =4. 

The Riemann - Hurwitz formula for the genus g of a n-sheeted covering of the Riemann 
sphere with the total multiplicity b of the branched locus gives 

b 

g = --n+l = 0. 
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A simpler way to establish rationality of the curve is to represent it as a quotient of a 
reducible rational curve 



[(s - l) 2 (s + 2) + 3w 2 ] [(a + l) 2 (s - 2) - 3u> 2 ] = 



over the involution 



(s, w) i ^ (— s, w). 
The projection of (1.3.32) to (2.1.19) reads 

(s, w) i— > (p = — 1 + s 2 , w). 

Further substitution 

s = 2 + 3t 2 

gives a rational parametrization of the second component 

(s + l) 2 (s-2) -3w 2 = 

of the curve (1.3.32) 

p = 3 + 12t 2 + 9t 4 ' 

> 

w = 3t(l + t 2 ). 
Now we can integrate the equation: 

, dw dt . . 

az = — = k, t = t&niz — zq) 

p 1 + t A 



that gives the solution in the form 



w = 3 tan 3 (z — zq) + 3 tan(z — zq) 



3sin(z — zq) 

COS 3 (z — Zq) 



Other components of the curve (1.3.32) give nothing new. 



(1.3.32) 



(1.3.33) 



(1.3.34) 



Exercise 1.3.8 Consider a differential equation of the form 

F(w',w,z) = 

where F(p, w, z) is a polynomial in 3 variables p, w, z. Assume that for any fixed g the genus 
of the algebraic curve F(p,w,z) = is equal to zero. Moreover, assume that the solutions to 
the differential equations have no movable critical singularities. Prove that the solutions to 
the differential equation expresses rationally through solutions to Riccati equation 



du 
dz 



a(z)u 2 + b(z)u + c(z) 



with rational coefficients a(z), b(z), c(z) 
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At the end of this section we consider an example of application of the above approach 
to a mult i- dimensional mechanical problem. 

The Kowalevskaya problem: consider the system of equations of motion of a rigid body 
with a fixed point 



Ax + {A 3 - A 2 ) uj 2 uj 3 = mg (c 2 7 3 - C372) 
A 2 L02 + (Ai - A 3 ) u 3 ui = mg(c 3 -fi - C173 

A 3 lu 3 + (A 2 - Ai)uiU2 = mg{c x ^2 - C271) , 

71 = w 3 72 - lo 2 73 

72 = ui 73 - uj 3 71 

73 = ^2 7i - ^i 72 J 



(1.3.35) 



(1.3.36) 



Here to = (ui, UJ2, oj 3 ) is the vector of angular velocity of the rigid body, 7 = (71,72,73) the 
coordinates of the unit vector along the Oz axis; the coordinates of both vectors are given 
with respect to the co-moving frame having the origin at the fixed point of the rigid body. 
The positive numbers Ai, A 2 , A 3 are the eigenvalues of the inertia tensor, c = (ci,C2,c 3 ) 
the (constant) coordinates of the baricenter, m is the mass of the body, g is a constant (the 
gravitational acceleration). What are the values of the parameters A\, A2, A 3 , a, c 2 , c 3 
for which the solutions (uj(t),^(t)) have no movable critical points in t G C? The following 
remarkable theorem was proved by S. Kowalevskaya (1889): 

Theorem 1.3.9 All singularities of solutions to the equations (1.3.35), (1.3.36) in the com- 
plex t-plane are poles if and only if 

(i) c\ = C2 = c 3 = (the case of Euler - Poinsot). 

(ii) A\ = A2 and c\ = C2 = (axial symmetry; the Lagrange - Poisson case). 
(Hi) A\ = A2 = A3 (spherical symmetry). 

(iv) A\ = A2 = 1A 3 , c 3 = (the Kowalevskaya case). 

In all four cases the dynamical equations (1.3.35), (1.3.36) reduce to a completely inte- 
grable Hamiltonian system. That means that, along with the Hamiltonian (the total energy) 



H = - (Aiul + ^2^2 + A 3 uj 3 ) - my (C171 + C272 + C373) 



(1.3.37) 



1 



M? Ml M\\ 

■> V AT + ~M + ^3 ' ~ mg ^ Cl71 + 0272 + C313 ' 



(here we introduce the vector of angular momentum 

M = (M 1 ,M 2 ,M 3 ), Mi = AiUi, i = l, 2, 3) 



(1.3.38) 
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and two geometric integrals 

(AT, 7) = Mi 71 + M 2 72 + M 3 73 

(1.3.39) 

(7, 7) = 7? + 72 + 73 

there is an additional independent first integral. Namely, in the Euler - Poinsot case the 
additional integral is 

I = (M, M) = Ml + M| + M|. (1.3.40) 
In the case of axial symmetry the additional integral is 

I = M 3 . (1.3.41) 

In the case of spherical symmetry the additional integral is 

(M, c) = c\M x + c 2 M 2 + C3M3. (1.3.42) 

The most nontrivial is the Kowalevskaya case where the additional first integral is quartic: 

/ = |A 3 (u;i +iu2) 2 +mg{c\ +ic 2 )(7i + * 72 ) | 2 • (1.3.43) 

Exercise 1.3.10 On the six- dimensional phase space with the coordinates (Mi, M 2 , M3, 71, 72, 73) 
introduce the following Poisson bracket: 

{Mi,Mj} = e ijk M k 

{M i , lj } = e ijklk (1.3.44) 
{7i,7j} = 

/ 1 2 3 \ 

Here eij k is the signature of the permutation ( ^ fc / summa ^ on over re P ea t e d index k 
is assumed. 

(i) Prove that the bracket (1.3.44) satisfies Jacobi identity. 

Hint: check first that linear functions on the phase space form a Lie algebra with respect to 
the bracket (1.3.44), namely, 

{(ai, M) + («i, 7 ), (02, M) + (a 2 , 7)} = («3, M) + (a 3 , 7) (1-3.45) 
03 = 01 x 02, a% = a\ x a 2 — a 2 x a\. 

Here the coefficients ai, ai of linear functions are considered as vectors inM 3 . The Lie algebra 
(1.3.45) is isomorphic to the semidirect product S03 x M 3 . 

(ii) Check that the functions (M, 7) and (7,7) are Casimirs of the Poisson bracket, i.e., 
that they commute with any other function on the phase space. 

(Hi) Prove that the equations of rotations of rigid body (1.3.35) is a Hamiltonian system 
with the Hamiltonian (1.3.37): 

M i = {H,M l }, j l = {H, 7l }, i = l, 2, 3. 
(iv) Check involutivity {H,L} = in all four cases (1.3.40) - (1.3.43). 
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According to the statements of this Exercise the level surface 

(M, 7 ) = fe, (7,7) = 1 (1-3-46) 

for an arbitrary constant k has a natural structure of symplectic manifolds. The restriction 
of the flow (1.3.35) to this level surface is a Hamiltonian system with two degrees of freedom. 
Thus it suffices to have one additional first integral to ensure complete integrability. 

In the classical cases (i) - (iii) the equations of motion (1.3.35) can be integrated in elliptic 
functions. Integration of the case (iv) performed by Kowalevskaya unraveled an important 
geometric structure behind this integrable case of equations of motion. Namely, consider the 
common level surface of the four first integrals 

(M, 7 ) =k 
(7,7) =1 

KS + ^ + 5)-(^ 1 + C2 7 2 ) =E 

|A 3 (oji +iuj2) 2 + mg(a + 102)^1 +ij2)\ 2 = K 

for generic values of the constants k,E,K. The dynamics of (1.3.35) takes place on this 
level surface. It turns out that a suitable compactification of this two-dimensional complex 
algebraic variety has a structure of an abelian variety, i.e., it is isomorphic to a 4-dimensional 
torus 

T 4 = C 2 /(integer lattice of rang 4). 

The dynamics becomes linear in the natural coordinates on this torus. The solutions to the 
Kowalevskaya problem can be expressed via theta-functions of two variables. The final form 
of these expressions was obtained by F.Kotter in 1893. 

The quite surprising connection of the property of absence of moving critical points with 
integrability of the system observed by Kowalevskaya was the starting point for the so-called 
Painleve test in the modern theory of integrable systems. 



(1.3.47) 



1.4 Poincare method of small parameter 

The method of small parameter, proposed by Poincare (1892) in his studies on celestial 
mechanics, is also an important tool in the study of properties of solutions of analytic differ- 
ential equations. The method is applicable to a system of differential equations depending 
on a parameter e, 

^=f(z,w,e) (1.4.1) 

where w = (w\, . . . ,w n ) is the vector of depending variables, the vector function f of the 
right hand sides is assumed to be analytic in z, w±, . . . , w n , e for sufficiently small |e|. The 
system (1.4.1) is considered as a perturbation of the so-called unperturbed system 

^ = f(z,w,0). (1.4.2) 
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Suppose we know how to solve the unperturbed system. What can be said about solutions 
of the perturbed one? 

Comparing with Cauchy theorem the novelty of the Poincare method is that, for systems 
depending analytically on a parameter in certain cases one can establish global existence of a 
solution. 

For the sake of simplicity we will formulate it for the case of a system of two first order 
differential equations, so we redenote w = (u,v), f = (/, g). 



Theorem 1.4.1 Consider a system of two differential equations of the form 



dv 
dz 



(1.4.3) 



g(z,u,v,e) 



with the right hand sides f(z,u,v,e), g(z,u,v,e) analytic in (z,u,v,e) for sufficiently small 
e. Assume that for e = the general solution 



u = (f>(z, Ci, C 2 ) 
v = il>(z,C 1 ,C 2 ) 





/ d<p> 


d<p 






' dd 


dd 


| , det 










dip 


dtp 


) 




\ dd 


dd 



7^0 



(1.4.4) 



to the system (1.4.3) is given. Moreover, assume that the solution (1.4.4) is analytic on a 
neighborhood of a curve 



C:z = z(t), t G [0,1], z(0) = z o , z(\) = z 1 
on the complex plane. Then there exists a unique solution 

u = (j)(z,Ci,C 2 ,e) 

> 

v = ^(z,C 1 ,C 2 ,e) j 



(1.4.5) 



(1.4.6) 



to the full system (1.4.3) analytic for sufficiently small e and also analytic in z on a neigh- 
borhood of the curve C . 



Proof: Step 1. Construct the solution as a formal series in e 

u = <f>-\- e<f>\ + e 2 4>2 + ... 

(1.4.7) 

v = V> + ei/ji + e 2 ip2 + 

We will prove that the coefficients are analytic in z near the curve C. They will be determined 
uniquely by the normalization 

cf)i(zo) = 0, i/ji(z ) = 0, i > 1. (1.4.8) 

Let us adopt the following short notations: denote /, g, df/du, df /dv, df/de etc. for the 
values of the functions f(z, u, v, e), g(z, u, v, e) and their derivatives at 

(z,u,v,e) = (z,4>,ip,0). 
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The construction: for the corrections (</>%, ipi) one derives from (1.4.3) a recurrent system 
of inhomogeneous linear differential equations of the form 



dz 



(1.4.9) 



where 



F k = F k 



d m f 



du % dvide l 



Gk = Gk 



d m g 



i+j + l= m<k, 
E.g., for we obtain the following system 



du l dv^de l 
a, b < k — 1. 



(1.4.10) 



#1 - 91 a , df I , df 
dz - du V 1 ^ dv VI ' de 



so F\ = df/de, G± = dg/de. For the second correction one has 



F 2 - 
G 2 



1 

2 

1 

2 



5 2 / , 5 2 /, 2 , df M , , d 2 f (2 



+ 



du 2 



L >{ + 2 



du dv 



(piipi + 



dv 2 



Pi 



+ ^4r<Pi 



d 2 g ^ d 2 g ^2^2 ^ 
de 2 du 2 1 dudv 



Mi 



d 2 g ... 



etc. Observe that the vector functions 
5ip 



and 



9^ 



du de 

d 2 g 
<9u <9e 



9C 2 

ac 2 



'1 + 



a 2 / 

<9t> <9e 
<9t> de 



(1.4.11) 



give a fundamental system of solutions to the associated homogeneous linear system 



Tz^= d /u^+"d^ 



9/ 



(1.4.12) 



due to the assumption (1.4.4). Therefore one can solve the system (1.4.9) iteratively by 
quadratures using the "variation of constants" method: 




dz' 



94>(z) d<t>{z) \ / 
dd dC 2 » 



ZD 



dip(z) dip(z) 



94>(z') 9<t>(z') \ 
dC\ dC 2 ' 



-1 



\ dd dC 2 / 




(1.4.13) 



The solutions are chosen in such a way that 

<Pk{z ) = t/Jk{z ) = 0, Vfe>l. 



(1.4.14) 
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By construction every term of the series is analytic on some neighborhood of the curve C. 
Step 2: proof of convergence. Denote 

F(z, 6<j>, Sip, e) := f(z, <P + S<P,iP + dip, e) - f(z, <j>, iP, 0) 

(1.4.15) 

G(z, Sep, Sip, e) := g(z, <p + S<p,ip + Sip, e) - g(z, p, ip, 0). 
By definition we have 

£Scf) = F(z,S<t>,SiP,e) 



d Sip = G{z, 5<p, 5ip,e). 



(1.4.16) 



dz 



Let 

F = Y J ^kim{z){S<P) k {SiP) l e m , a 000 (z) = 

G = hi m (z)(S<P) k (Sip) l e m , b oo(z) = 

be the Taylor expansions of the analytic functions F(z,S(p,5ip,e), G(z,5(p,5ip,e). Choose a 
small positive number p in such a way that 

\F(z,S(f>,Sip,e)\ < M ) for z G C 

\sp\ < P 
\Sip\ < p 

\G(z,S<p,Sip,e)\ <M J |e| < p. 
From the Cauchy inequalities 

t M M 

\a k im{z)\ < pk+l+m 

il r M M 
\hlm[Z)\ < 



pk+l+rn 



one readily derives the following majorants for the functions F and G: 
F(z, Sep, Sip, e)±Yl (ScP) k (6iP) l e m - M 



M ty + fy + e ^ M (S^ + SiP + e) (l + #±g±^) 

P P 



same for G. 

Consider an auxiliary system 



d M W + ^ + £) (i + *±J**) 

— 00 = 

dz p _ <5<^+<5V'+<; 



P 



(1.4.17) 



d_ M (^ + # + e )h+ /; 

— OW = 

P 
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Let us find the solution to this system with the initial data 

S(f)(zo) = Sip(zo) = identically in e. 
Due to the symmetry of the system and of the initial condition with respect to the permutation 

5(f) <-> 5if> 

for the solution the identity 



holds true. Denote 



A := 



5(f){z) = 6if>(z) 
5(f> + 5i) + e 25~4> + e 



P P 

For the function A(z) one obtains the following Cauchy problem 

dA M 1 + A a / \ e 

-f- = 2 — A -, A(zq ) = -. 

dz p 1 — A p 



(1.4.18) 



Separating the variables 



2 M^ = (l-A) rf A 
p A(l + A) 



^A 1 + A 

and integrating one obtains the solution in the form 



dA 



A 



,-e p 



(z-zo) 



(1 + A) 2 (e + pY 
Denote a the right hand side of the equation. For z = zq 



a(z ) 



(e + p) 



2 4 



Solving the quadratic equation 



A = -1 + — [1 ± y/l - 4aJ 

Zi(X 



yields 



A(^ ) 



■1 



(e + p) ; 
2ep 



1 ± 



P~ e 
P + e 



, for + sign 
, for — sign. 



Therefore one must choose the "-" sign. Finally one obtains the solution in the form 

(e + p) 2 



A(z) 



-1 + 



2ep 



-e p ' 



2^2 



1 



1 



ep ^(z-z ) , o e "P 



7 e p 



(e + p) 2 
This function is analytic for 



+ 2 



(e + p) 



4ep gM (z _ Zo) 
(e + p) 2 



(e + p) f 



4ep c Mi ( ,_ 2o) 



< 1, 



(e + p) 2 

i.e., for sufficiently small |e| for any z near the compact C. Like in Section 1.1 from analyticity 
of solutions to the majorant system (1.4.17) one deduces convergence of the "perturbation 
series" solutions (1.4.7), (1.4.8) to the initial system. □ 
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1.5 Method of small parameter applied to Painleve classification problem 

Painleve addressed the problem of classification of the second order ordinary differential 
equations 

w" = R(w',w,z) (1.5.1) 

having no movable critical singularities. His application of the Poincare method to this 
problem is based on the following simple corollary from the Theorem 1.4.1. 



Theorem 1.5.1 Under the assumptions of Theorem 1.4-1 consider the case of the closed 
curve C, z\ = zq. The solution (u(z, e), v(z, e)) is single-valued along C iff the coefficients 
<f>k(z), ipk( z ) of the perturbative expansion (1.4.7) are single-valued for all k>0. 

Proof: Let the functions (f>k(z), ipk( z ) satisfy 

(pk(zi) = <Pk{z ), tpk(zi) = ipk(zo) for k < K, 



f>K(zi) - (Pk(z )\ 2 + \iPk(zi) ~ iPk(z )\ 2 / 0. 



Then 



lim e 

€^0 



-2K 



\4>(zi,e) - (f>(zo,e)\ +\ip(z!,€) -i/>(zo,e)\ / 0. 



□ 



To warm up let us explain how to obtain a simple proof of Theorem 1.2.13 by applying 
the small parameter method. We will prove first that the rhs of the equation 

w = f(w,z) 

must have no poles. In the opposite case one can locally represent the equation in the form 

w'= , F{W \\ k , F(a(z),z)?0 
(w — a(z)) K 

with a positive integer k. After the substitution 

w = a(z) + u 

the equation will read as 

, _ G(u,z) 

u — u 
U K 

for some function G(u, z) such that G(0, zq) / for some zq. Observe that after the substi- 
tution 

z = z + e m T 



u = e v 



the equation will take the following form 



du = ^ (k+1)m G(e n v, z + e m r) 
dr u k 
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The choice 

m = 1, n = k + 1 

yields the following equation 

du = G(e n v,z + e m T) = G(0,z ) + 
cZt u fc u fc 

with the rhs depending analytically on the small parameter e. The solutions of the unper- 
turbed equation 

du _ G(0, zq) 
dr u k 
have critical singularity at z = zq if k > 0: 

i i 
u=8thi, a = [(k+ l)G(0,z )] fc+1 • 

Since zq can locally take arbitrary values, the singularity will be movable. So, to avoid critical 
singularities one must have k = 0. 

The next step is to prove that, for a polynomial rhs the equation 

w' = ao(z)w h H h afc(z) 

has no critical singularities only if k < 2. The susbtitution 

z = z + e m r 
tt; = e~ n u 



gives 

du 
Ik 

The choice 



^m-n(fc-l) 



a (zo + e"V) n fc + e n ai (z + e m r) n^ 1 + • • • + e nk a k (z + e m r) 



n = l, m = k — 1 
yields the following unperturbed equation 

— = a (z )u . 

The solution 

u = — b = [(1 - fc)a (zo)] ^ 

7- fc-1 

has a critical singularity if k > 2. 

Let us return to the Painleve classification problem. We are looking for the second order 
ordinary differential equations of the form 

w" = R(w',w,z) (1.5.2) 

admitting no movable critical singularities (the so-called Painleve property). It is assumed 
that the right hand side R(w', w, z) is a rational function in (w', w) with coefficients that are 
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analytic function in z on some domain 0, C C. There are some obvious classes of equations 
satisfying Painleve property. First, second order linear differential equations with rational 
coefficients are already known to meet the Painleve condition. Next, second order differential 
equations of the form 



w" = 3a(z) w w' 



p(z) + 2 



/ x / \ a '{ z ) a"(z) 



a'(*) 
a{z) _ 



w + 



w' - a 2 (z)w 3 + [a(z)p(z) + 3a'(z)] w 2 



r(z) 



(1.5.3) 



a(z) 



a(z) a(z) 
obtained from the linear equation 

y"' + p{z)y" + q{z)y' + r(z)y = 
with rational coefficients by substitution 

i y' 



w = — 



a(z) y 



(cf. (1.2.6)). One has to also take into account second order differential equations solved by 
elliptic functions. E.g., differentiating Weierstrass equation 

{p'f =4p 3 -92P-92, 
in z and dividing by 2p' one arrives at 



P =®P - 7T- 



(1.5.4) 



Remarkably, besides these essentially obvious classes of second order differential equa- 
tions satisfying Painleve property there are six new equations discovered by Painleve and his 
student B.Gambier. They are the celebrated Painleve equations: 



w = 



w 



w = 



w 



w 



w 



6w 2 + z 

+ z w + a 

a w 2 +f3 



UJ W I 

W Z 



II _ w' 



2it] + | w 3 + Az w 2 + 2{z 2 - a)w + 



2w ~ w-1 



l2 w < . (w-1) 2 / . 3 , 



■yw _|_ Sw(w+l) 
z ' w—1 



. w(w— l)(w— z) 

2 2 (Z-1) 2 



n ,3z, 7 (z-l) , Sz(z-l) 
wj 2 " (w-1) 2 (w-z) 2 



Pi 
Pii 

Pm 

Piv 

Pv 

Pvi 
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Theorem 1.5.2 (i) Solutions to the differential equations Pi - Pyi have no movable critical 
singularities. 

(ii) General solutions to these equations are new transcendental functions, i.e., they can not 
be expressed via algebraic functions, or via solutions to linear differential equations, or via 
elliptic functions. 

(Hi) Given a second order differential equation of the form (1.5.2) with the right hand side 
depending rationally on w, w' with coefficients analytic in z, suppose that solutions to this 
equation satisfies Painleve property. Then solutions to this equation can be expressed via 
algebraic functions, or elliptic functions, or solutions to linear differential equations, or via 
solutions to one of the equations Pi - Pyi- 



In this section we will give some hints about the proof of the last statement of the theory, 
namely, we develop an approach to the classification of the second order differential equations 
satisfying Painleve property based on the small parameter method. 

Proof: Step 1. We prove that the rhs R(w',w,z) must be a polynomial in w' of degree at 
most 2. 

Let us first prove that poles in w' are not allowed in the rhs. Indeed, in the opposite case 
the equation near the pole reads 

n = Q(w',w,z) 



[w' — a(w, z)] k 

where k is a positive integer and Q(a(w, z),w, z) is not an identical zero. The substitution 

u = w' — a(w, z) (1.5.5) 

represents the equation as a system 

dw 

— — = a(w, z) + u 
dz 

du Qi(u,w,z) 
dz u k 

with some new function Qi(u, w, z) analytic and non vanishing near u = 0. Introducing a 
small parameter 

z = z + e p Z, w = w + e q W, u = e r U 



with arbitrary (zq,wo) yields 



^-pdK = a ( WQ + e i W: ZQ + € p Z ) + e r U 

ClZi 

(k+i)r-pW = Qi{e r U,w + eiW,z + ePZ) 

dZ u k 



The choice 

q = p = (k + l)r 
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of positive numbers p, q, r gives the following unperturbed system 

dW 

— = a{w ,z ) 
dU _ Qi(0,wo,z ) 

dZ ~ u k 

Integration of the last equation gives critical points for U(Z) if the integer k is positive: 

U =[(k+l)Q 1 (0,w o ,z o )Z + C]kh . 



However, from the definition (1.5.5) it follows that the function u{z) cannot have critical 
points if the solution w(z) has none. 

We proved that the rhs of (1.5.2) is a polynomial in w'. Let us show that the degree of 
this polynomial is at most equal to 2. Indeed, rewriting the differential equation 



w" = Aq(w, z)w' n + Ai(w, z)w' n 1 + • • • + A n (w, z) 



as a system 



dw 
dz 

^ = A (w, z)u n + At (to, z)u n ~ x H V A n (w, z) 

introduce the small parameter as follows 

z = z + e p Z, w = w + e q W, u = e~ r U. 

The system will read 

q+r-p _ 

dZ 

e r(n-l)-p dR = + ^ Z() + £ p ^jjn + 0(e) 

dZ 

Choosing 

p = q + r = r(n — 1) 
gives the following solution to the unperturbed system 

U = [C - (n - l)A (w ,z )Z}^ . 

This is a critical pole if n > 3. We arrive at an equation of the form 

w" = A {w, z) w' 2 + Ai(w, z) w' + A 2 (w, z) (1.5.6) 

with some rational functions Aq(w,z), A\(w,z), A2(w,z). 

In these notes we will be unable to reproduce all details of the proof of the Part (iii) of 
the Painleve - Gambier theorem. However, to give some flavour of the remaining part of the 
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proof let us jump to considering the particular case of equations of the form (1.5.6) having 
Aq = A\ = 0. So, consider an equation of the form 



w" = R(w, z) 

satisfying the Painleve property. As before, the function R(w, z) is rational in w. 
Step N, (N >> 1). Let us prove that all poles of the rhs in w are simple. 
We will now proceed to studying the particular case of equations of the form 

w" = ao(z)w 3 + a\{z)w 2 + a2(z)w + a^{z). 

First of all, in the case ao(z) / doing the transformation 

z >-> z = f(z) 



(1.5.7) 



(1.5.8) 



(1.5.9) 



w ^ w = A 1 w 



with 

one can achieve 



f'(z) = al /3 (z), X(z) = V2a Hz) 



a (z) = 2. 

If this is already the case in (1.5.8) then, after the shift 



w ^ w — - ai(z) 



the equation reduces to the form 



w " = 2w 3 + a(z)w + b(z). 



(1.5.10) 



We have now to derive conditions for the functions a(z), b(z) necessary for absence of movable 
critical singularities. 

Introducing a small parameter by means of the substitution 

w = <r x W, z = ZQ + eZ 

we obtain 

W" = 2 W 3 + e 2 a W + e 3 [a' Z W + b ] + e 4 
Here we use the following short notations 



a'^Z 2 W + b' Z 



+ 0(e 5 ). (1.5.11) 



The unperturbed equation 
has a first integral 



a :=a(z ), a := a'(z)\ z=zo 

W" = 2W 3 
W' 2 = W 4 + C. 



etc. 
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The general solution to this equation is given in elliptic functions. However, for the particular 
case C = it has a simple solution 

Now let us look for the perturbative series solution 

W = + e 2 2 + e 3 (/>3 + e 4 04 + .... (1.5.12) 
After the substitution one obtains the following recursive equations 

02 - 6 (f> 2 4>2 = a <f> 

<p3-6<j) 2 (p3 = a Z ( j ) + b (1.5.13) 
4>'[ - 6 <f> 2 (p 4 = X - a' ' Z 2 $ + b' Z + 6(jxf> 2 2 + a (j)2- 

Choose first 



T Z + C 

The following statement will be useful in sequel. 

Exercise 1.5.3 Linear differential equation of the form 

x n y (n) + aix n ~ l y( n ~ l) + ■■■ + a n -ix y' + a n y = 0, ai, . . . , a n _i, a n G C (1.5.14) 
is called Euler equation. The polynomial 

P(k) := k(k - 1) . . . (k - n + 1) + a ± k(k - 1) . . . (k - n + 2) H h a n _i/c + a n (1.5.15) 

is called the characteristic polynomial o/ t/ie Euler equation. 

(i) Given a root ko of the characteristic polynomial, 

P{k ) = 0, P'(k ) + (1.5.16) 

prove that 

y(x)=x ko (1.5.17) 

is a solution to the Euler equation. 

(ii) If ko is a root of multiplicity m, i.e. 

P(k ) = 0, P'(k ) = 0,...,P^ 1 \k )=0, P^/O (1.5.18) 
then prove that the functions 

yi = x h \ y 2 = x k Hogx,...,y m = x kQ log m - 1 x (1.5.19) 

are linearly independent solutions to the Euler equations. 

(Hi) Prove that solutions of the form (1.5.17), (1.5.19) where ko runs through the set of 
all roots of the characteristic polynomial give a basis in the space of solutions to the Euler 
equation. 
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The linearized equation 

(Z + C)V-6Y> = (1.5.20) 
is a particular case of Euler equation. The characteristic equation 

fc(A; — 1) — 6 = 

has two roots 

k = -2 and k = 3. 
Thus the general solution depending on two arbitrary constants k±, k 2 reads 



^ = {z + cf +k2{z + c) ■ (L5 - 21) 

The following statement will be of use for solving inhomogeneous Euler equation. 



Exercise 1.5.4 Consider inhomogeneous Euler equation of the form 

xV n) + aix™- 1 ^™- 1 ) + • • • + a n . lX y' + a n y = ax K . (1.5.22) 

(i) Prove that, if k is not a root of the characteristic polynomial (see Exercise 1.5.3 above), 
P(k) 7^ 0, then a particular solution to the inhomogeneous equation is given by the formula 

" x K . (1.5.23) 



" P(K 

(ii) If k is a simple root of the characteristic polynomial, i.e. 

P(k) = 0, P'(k) + 
then prove the following formula for a particular solution: 

y=——x K \ogx. (1.5.24) 

f [K) 

Generalize to the case of roots of higher multiplicities. 

Applying this method we find general solutions to the linearized equations 

= ~\ a {Z + C) + {z ^ c)2 +k 2 {Z + Cf (1.5.25) 

03 = I a' C(Z + C)-\ K + b ] (Z + Cf + {z ^ c)2 +k 2 (Z + Cf. (1.5.26) 
Integration of the subsequent equation yields logarithmic terms: 

04 = ^ Z +Cf Ko+26'o] log[Z + C]-^ K' + 26 / ] {Z + Cf + \ C « + &' ) - l -C 2 a'^Z+C) 

(1.5.27) 

(we choose integration constants k\ = k 2 = in (1.5.25)). Cancellation of logarithmic terms 
implies 

a"(zo) + 2b'(z ) = 0. (1.5.28) 
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Doing similar calculations with the perturbative solution starting with 
yields another constraint: 

a"(z ) - 2b'(z ) = 0. (1.5.29) 

Therefore 

a"(zo) = 0, b'(zo) = 0. 
Since zq is an arbitrary complex number we finally obtain that 

a(z) = Xz + n, b(z) = v (1.5.30) 

with some constants A, fi v. For A = one obtains differential equation solved in elliptic 
functions. If A 7^ then, doing if necessary an affine transformation of the independent 
variable z and a rescaling of w we finally arrive at the Painleve-II equation 

w = 2w + zw + a. 

In a similar way one can deal with equations of the form 

w" = ao(z)w 2 + ai(z)w + (i2{z), ao(z) / 0. (1.5.31) 
After the substitution of the form (1.5.9) with 

\ = 6a Hz), f{z)=al l \z) 
and a suitable shift of the dependent variable one reduces the equation to the form 

w" = 6w 2 + a(z). (1.5.32) 
Introducing the small parameter by the substitution 

z = z + eZ, w = e~ 2 W 
one arrives at the following equation 

= 6 W 2 + £ 4 ao + £ 5 a / z + 1 e 6 a // Z 2 + Q ^ (1 5 33) 

The notations ao, Oq etc. are similar to those used above. 

For the unperturbed equation choose the particular solution 

m = TzTcf 

The coefficients of the perturbative solution 

W = 4>(Z) + e 4 MZ) + e 5 MZ) + ^M^) + ■■■ 
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are determined from the following recursive procedure: 

04 - 12 04 = a 
5 '- 12005 = a' Z 

06 - 12 006 = o a'o'Z 2 . 



The general solution of the homogeneous linearized equation 

(Z + C)V -12^ = 

is obtained in the form 

Applying the method of Exercise 1.5.4 we obtain single valued solutions 04 and 05 but 



b & = ±-<%{Z + Cf log(Z + C) - \alC (Z + C)- lag + h(Z + C) 4 + 
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14" uv ' ~> m ' ~> 3 u v ' ; 12 " u ' " v ' ~ ; ' {Z + Cf 
This function has no critical point at Z = —C iff 

a"(z ) = 0. 

Due to the freedom in the choice of zq the function a(z) must be linear, 

a(z) = Xz + /x. 

In the case A = the differential equation (1.5.32) integrates in elliptic functions. Otherwise 
doing an affine change of the variable z and a rescaling of w (1.5.32) becomes the Painleve-I 
equation 

w" = 6 w 2 + z. 

□ 

It remains to prove 1 that solutions to the equations Pi - Pyi have no movable critical 
singularities. Actaully we will explain in sequel how to prove that 

• solutions to Pi, Pn and Piy equations are meromorphic functions on C having essential 
singularity at the fixed singularity z = oo; 

• solutions to Pin and Py are meromorphic functions on the universal covering of C* = 
C\{z = 0}; 

• solutions to Pyi are meromorphic functions on the universal covering of Riemann sphere 
with three punctures (the fixed singularities) z = 0, z = 1 and z = oo. 

In all the cases the position of fixed singularities can be easily guessed by looking at the 
singularities of coefficients of the equation. The proof of meromorphicity is more involved. 



1 One has to also prove that general solutions to the equations Pi - Pvi cannot be expressed via known 
functions. The proof of this statement required some technique of differential Galois theory. It will not be 
given in the course. Besides, we have to note that some particular solutions to these equations can be expressed 
via known special functions, some of them even via elementary functions. The classification problem of these 
particular solutions can be difficult. For example the problem remains open for the case of Painleve-VI 
equation. 
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The technique to be used is based on a realization of the Painleve equations as monodromy 
preserving deformations of linear differential equations with rational coefficients. In the next 
section we will begin with explaining the basics of the monodromy theory of solutions to 
these equations. 
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2 Part 2: Monodromy of linear differential operators with 
rational coefficients 



2.1 Solutions to linear systems with rational coefficients: the local theory 



Let us begin with recalling some basics from the theory of systems of linear differential 
equations. Consider an n-th order system linear homogeneous system of the form 



y' = A( z ) V: y = 



( vi \ 

y2 



\y n J 



( anO) 

021(2) 



, A(z) = 



ain(z) \ 

(I2n{z) 



, zeC. 



(2.1.1) 



\ a nl (z) . . . a nn {z) J 



The space of solutions defined in a neighborhood of a regular point zq of coefficients aij(z) is 
a linear space of dimension n. Indeed, there exists a unique solution to (2.1.1) with arbitrary 
initial data 

/ y°i \ 



yi(zo) = yi, y2{z ) = 2/2, • • • ,yn(zo) = y°, 



y° 2 



\y°nJ 



g c n . 



Choosing a basis of n linearly independent solutions one can represent an arbitrary solution 
as a linear combination 



y{z) = a 



( Vu(z) \ 

V2l{z) 



\ Vnl(z) J 



H h c n 



( yin(z) \ 

y2n(z) 



\ ynn(z) J 



(2.1.2) 



with some constant coefficients ci, . . . , c n . The basic vector functions are columns of the 
fundamental matrix for the system (2.1.1) 



Y(z) 



( yu(z) ■ ■ ■ yin(z) \ 



\ Vnl(z) . . . y nn {z) J 



Linear independence of the vector functions is equivalent to never vanishing of the determi- 
nant of this matrix 

W(z) :=dety(z) 7^0. (2.1.3) 
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Any solution (2.1.2) can be represented in a matrix form as follows 



C2 



y(z)=Y(z)c, c 



(2.1.4) 



V Cn J 

with an arbitrary constant vector c £ C ra . 

The fundamental matrix is a matrix solution to the system (2.1.1): 

Y' = A{z) Y 



(2.1.5) 



Conversely, any matrix solution to (2.1.5) satisfying (2.1.3) is a fundamental matrix of the 
system. 

Lemma 2.1.1 Any two fundamental matrices Y{z), Y{z) for the system (2.1.5) are related 
by a linear transformation of the form 



Y{z) = Y(z)C 
with a nondegenerate constant n x n matrix C . 

Proof: follows from (2.1.4). 



(2.1.6) 



□ 



Example 1. Linear systems with constant coefficients 

y' = Ay. (2.1.7) 

A fundamental matrix is given by the matrix exponential function 

Y(z) = e Az . (2.1.8) 

Let ai, . . . , a n be the eigenvalues of the matrix A (more precisely, the roots of the charac- 
teristic polynomial of A). Denote 

a := diag(ai, ...,a n ) 

the diagonal matrix. Suppose that the matrix A is diagonalizable. For example this is always 
the case if the roots of the characteristic polynomial are pairwise distinct. Denote T the 
matrix of eigenvectors of A. The conjugation by this matrix reduces A to the diagonal form 

T~ l AT = a. 



Then a fundamental matrix can be computed as follows 

y(z) = Te", where e &z = diag (e aiZ , . . . , e a " 2 ) . 



(2.1.9) 
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It is related to (2.1.30) by the multiplication by T 1 on the right. In the presence of multiple 
roots of the characteristic polynomial 

det (A - a ■ Id) = 

the fundamental matrix (2.1.30) involves also polynomials. They can be computed by reduc- 
ing the matrix A to Jordan normal form and using the following formula for exponential of 
an elementary Jordan block of size k: 



e Jz = e az 



( 1 







V o 



2! 



3! 
2! 



1 





(fc-2)! 



Z 
1 



J 



(a 






V 



1 

a 




1 








\ 





a 1 
a J 

(2.1.10) 



Example 2. The system 



y = —y 

z 



(2.1.11) 



with constant matrix A can be reduced to (2.1.7) by the substitution 

1 dy 

logz = t, ^ = Ay. 

Thus, for a diagonalizable matrix A the fundamental matrix (2.1.32) becomes 

Y = Tz & , where z & = diag (z a \ . . . , z an ) , (2.1.12) 
01, . . . , a n are eigenvalues of the matrix A, T~ x AT = a. Observe that the function 



is a multivalued function on C* = C \ {0} unless \i is an integer. However, the following 
simple statement holds true. 

Lemma 2.1.2 For an arbitrary fj, G C there exists an integer N and a positive p such that 
for an arbitrary choice of the branch of the function z^ and arbitrary two real numbers a < (5 
there exists a constant C such that 



z^\<C\z\ N for \z\ < p, a<aigz<(3. 



(2.1.13) 



Proof: Denote /Ur, fj,\ respectively the real and imaginary parts of p. Then 



\ Z P\ — e (/*R+»w)(log|«|+i argz) _ argzi 
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We choose 



I = max{/c £ Z \ k < p^} 
= 1 

f e-««, mi > 
" \ e~^, tn < 



iV = 



□ 



To be short we will say that the multivalued function has a polynomial growth at the 
singular point z = 0. It is important to impose restrictions for the argument: if z — > but 
is not constrained within a sector on the complex plane then the polynomial growth can be 
violated. 

Like in the previous example for a nondiagonalizable matrix A logarithmic terms will 
appear in the fundamental matrix (cf. the theory of Euler equation above). The solutions 
will still have a polynomial growth at the origin. 

Example 3. Consider now equation 

y' = ^y, rez, r > 2. (2.1.14) 



This case can also be reduced to a system with constant coefficients by a substitution 

* = -^7^I, § = Ay. (2-1.15) 

r — 1 z r 1 at 

So the fundamental matrix, for a diagonalizable matrix A can be chosen in the form 

Y(z) = Te~^ zl ~ r (2.1.16) 



in the notations of the previous examples. The matrix entries of this fundamental matrix are 
single valued functions on C* but they have an essential singularity at z = 0. We leave as an 
exercise to establish the structure of the fundamental matrix in the case of a nondiagonalizable 
matrix A. 



We will now address the problem of studying the local behaviour of solutions to linear 
systems of differential equations near a singular point of coefficients analytic on a punctured 
disk. To be more specific we will assume that the coefficients have a pole at the puncture. 
Without loss of generality we will assume the puncture to be put at the origin. The system 
then will be written in the form 

z r y' = A(z)y (2.1.17) 
for some integer r > 1 where the matrix valued function A(z) is analytic on the disc 

\z\ < p 

for some positive p and 

A(0) ± 0. 
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Definition 2.1.3 The number r — 1 is called the Poincare rank of the singular point z = 0. 
The singularity of the system (2.1.17) will be called Fuchsian if r = 1. 

We have seen, looking at simple examples, that the local behaviour of solutions near 
the singularity is qualitatively different for the cases r = 1 (at most polynomial growth in 
the sense of Lemma 2.1.2) and r > 1 (exponential growth). We will now analyze the local 
behaviour for the more general case of systems (2.1.17). We will see that the treatment is 
different for Fuchsian and non- Fuchsian cases. One can see some difference by applying the 
method of small parameter to the systems 

zy' = (A + A 1 z + A 2 z 2 + ...)y (2.1.18) 

or 

z r y' = (A + A lZ + A 2 z 2 + ...)y, r > 2. (2.1.19) 
Introducing the small parameter by rescaling 

z i ^ ez 

one obtains a regular dependence on e in (2.1.18) 

z y = (A + e A\z + e 2 A 2 z 2 + . . . ) y 
but a singular one in (2.1.19) 

z r y' = — [ (A + eA 1 z + e 2 A 2 z 2 + ...)y. 

An equivalent approach to the study of singularities of the fundamental matrix is based 
on the theory of gauge transformations of the system (2.1.17). 

Given an invertible matrix valued function G(z) defined on the punctured disk, one can 
transform the system (2.1.1) by the substitution 

y(z) = G{z)y{z) (2.1.20) 

to another system of the same type 

y' = A(z)y (2.1.21) 

with 

A = G~ l AG- G~ l G'. (2.1.22) 
The transformations of the form (2.1.22) are called gauge transformations. 

Definition 2.1.4 Two systems of the form (2.1.1), (2.1.21) are called strongly equivalent if 
there exists a gauge transformation (2.1.22) with the matrix valued function G(z) analytic on 
some disk \ z\ < p with never vanishing determinant. They are called weakly equivalent if such 
a gauge transformation with an analytic matrix valued function exists with the determinant 
never vanishing for \z\ > 0. 

The proof of the following statement is straightforward. 
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Lemma 2.1.5 Given a system of the form (2.1.1), then any weakly equivalent system has a 
form 

z f y' = A(z)y, A(z) + (2.1.23) 

with some r > 0, or the new system has no singularity at z = 0. Any strongly equivalent 
system has the form (2.1.23) with r = r. 

Let us begin with the Fuchsian case. 

Definition 2.1.6 A square matrix M is called resonant if there exist two eigenvalues Ai, A2 
of M such that the difference X\ — A2 is a positive integer. If such a pair does not exist then 
the matrix M is called nonresonant . 

The main result of the local theory of Fuchsian systems is the following 

Theorem 2.1.7 Any Fuchsian system of the form (2.1.18) with nonresonant leading term 
Ao is strongly equivalent to the "unperturbed" system 

zy' = A y. (2.1.24) 



Proof: We have to find an invertible analytic matrix valued function G = G(z) satisfying 
the following linear differential equation 

zG' = A{z)G - GA . (2.1.25) 

The proof will consist of two steps. First we will construct a formal series solution 

G(z) = G + zG 1 + z 2 G 2 + ... (2.1.26) 

to (2.1.25). At the second step we will prove convergence of the series. 

For the sake of technical simplicity the proof of the first part will be done under the 
additional assumption of diagonalizability of the matrix Aq. Doing if necessary a conjugation 

A(z) ^T- 1 A{z)T 

with a suitable invertible constant matrix T we may assume the matrix Aq to be diagonal, 

A = diag(Ai, . . . , A n ) =: A. 

Then we choose Go = Id. Substituting (2.1.26) to (2.1.25) and collecting the coefficient of z k 
we obtain a recursion relation for the coefficients G±, G2, 

fc-i 

[k,G k ]-kG k = -A k -Y J A i G k _ i , k = l,2,... (2.1.27) 

i=l 

Assuming the coefficients G±, . . . , G k -i already known we have to determine the matrix G k - 
Observe that for a n x n matrix M = (M^) the commutator [A, M\ = AM — MA has the 
following entries 

([A,M]).. = (A i -A j )M ij 
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So we can put 

<«*>« = -^d^(^ + |> G *--)- 

Due to the nonresonancy assumption the denominators Aj — Xj — k never vanish. 

We have proved, under assumptions of nonresonancy and diagonalizability of Aq, existence 
of a formal gauge transformation between the systems (2.1.18) and (2.1.24). The assumption 
of diagonalizability is purely technical and can be eliminated due to the following statement. 

Exercise 2.1.8 Let A be an arbitrary n x n matrix. Consider the linear map 

Mat(n, C) -»• Mat(n, C), M ^ [A, M] (2.1.28) 

of the space C™ 2 = Mat(n, C) of all n x n matrices to itself. Prove that characteristic roots 
of this map belong to the set of all differences Aj — Aj where X\, . . . , X n are the characteristic 
roots of the matrix A. 

We now proceed to the proof of convergence of the series (2.1.26). This will follow from 
the following general statement. 

Lemma 2.1.9 Consider a system of N linear differential equations 

z y = A(z)y 

with a Fuchsian singularity at the origin and the right hand side analytic for \z 
the system has a formal series solution 

y( z ) = y + z y 1 + z 2 y 2 + .... 
Then the series converges for \z\ < r. 

Proof: Denote Ao, Ai, A2, ■ ■ ■ the coefficients of Taylor expansion of the analytic function 
A(z), 

A{z) =A + zAi + z 2 A 2 + ■■■■ (2.1.31) 
For the coefficients of the formal series solution to (2.1.29) one obtains the following recursion 

k 

A y = o, ky k = A Q y k + Y J A i y k - l , k>i. (2.1.32) 

1=1 

Because of convergence of the series (2.1.31) there exist positive constants c and p such that 

\\A k \\<cp k , k = l,2,... 

that is, the inequality 

\A k y\<cp k \y\ (2.1.33) 
holds true for an arbitrary vector y. Indeed, it suffices to choose 

1 

p>-. 
r 



(2.1.29) 
< r. Suppose 

(2.1.30) 
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Choose a sufficiently big positive integer m in such a way that 

\A y\<m\y\ 

for any vector y. It is easy to see that the matrix Aq — Kid is invertible for any k > m since 

\A y - ny\ > \k \y\ - \A y\\ >(n-m)\y\^o vy^o. 

Thus the linear inhomogeneous equation 

A y-ny = z 

with an arbitrary right hand side Z for any k > m has a solution. Moreover for the norm of 
the solution one has an estimate 

\y\ < 

k — m 

A somewhat stronger estimate 

\y\ < \z\ 

follows if k > m + 1 . 

We want to establish estimates for the coefficients determined by (2.1.32) 

|34|<Ci? fc (2.1.34) 

for some positive constants C, R. More precisely, we first choose sufficiently large R > p in 
such a way that 

Cf> <1. (2.1.35) 



R-p 

Next, choose the constant C in such a way that the inequalities (2.1.34) hold true for small 
values of k, namely, for < k < m. 

Let us now derive (2.1.34) inductively for all k. Rewriting (2.1.32) in the form 

k 

Aoyk - ky k = z k , z k = - ^2 Aiy k -i 

i=i 

and using (2.1.33) and the inductive assumption we obtain 

k k—1 . oo 

\Z k \ < Y,cp*CR*-> = cCpR^ £ (g < cCpR^Q' = CR k J^~ p 

i=l i=0 i=0 ' 

Due to the choice (2.1.35) we derive that 



\Z k \ < CR k . 

Finally, for k > m + 1 this implies that \y k \ < C R k . 
We proved that the series (2.1.30) converges for 

1 
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The last step in the proof of Lemma is in observing that any analytic solution to the lin- 
ear system (2.1.29) can be analytically extended onto the entire domain of analyticity of 
coefficients, i.e., the series (2.1.30) converges for all \z\ < r. □ 

We can now complete the proof of the Theorem. We can consider the equation (2.1.25) as 
a system of linear differential equations for the n 2 -dimensional vector valued function G(z). 
We have proved above existence of a formal series solution to this system. Due to Lemma the 
formal series converges. Since G(0) = Go = Id the matrix G{z) is invertible for sufficiently 
small \z\. □ 



Corollary 2.1.10 A Fuchsian system (2.1.18) with diagonalizable nonresonant matrix Ao = 
A(0), 

T- 1 A)T = A = diag(Ai,...,A n ) 
has a fundamental matrix of the form 

Y{z) = G{z)z A , z A = diag (V 1 , . . . , z A ") (2.1.36) 

with the matrix G(z) analytic and invertible for sufficiently small \z\. 

What happens in presence of resonances? We give an answer leaving the proofs as an 
exercise. 

Let J be the Jordan normal form of the matrix Aq. Recall that it is an upper triangular 
matrix determined uniquely up to a permutation of Jordan blocks. We will choose a good 
Jordan form in such a way that the diagonal entries Ai, . . . , X n (not necessarily distinct) of 
the matrix J satisfy 

Re Ai > Re Xj for i < j. (2.1.37) 

Decompose 

J = A + R (2.1.38) 

where 

A = diag(Ai,...,A„) 

is a diagonal matrix and Ro is a collection of Jordan blocks with all zeroes on the diagonal. 
The upper triangular matrix Rq is nilpotent. The only non-zero entries of this matrices are 

(R )..^0 only if i^j and Aj = \j. (2.1.39) 

Theorem 2.1.11 There exist upper triangular matrices R±, R2, ■ ■ ■ such that 

{Rk^j^O only if X i -X j = k (2.1.40) 

such that the Fuchsian system (2.1.18) is strongly equivalent to a system with polynomial 
coefficients of the form 

zyf = (A + R + zR 1 + z 2 R 2 + ...)y (2.1.41) 
(only finite number of non-zero terms). 
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Theorem 2.1.12 The system (2.1.41) has a fundamental matrix of the form 

Y(z) = z A z R , R = R + R 1 + .... (2.1.42) 

Hint: Use the following property of the matrices R k '- 

z A R k z~ A = z k R k , k = 0, 1, 2,.... (2.1.43) 



Corollary 2.1.13 Any system of linear differential equations near a Fuchsian singularity 
has a fundamental matrix of the form 

Y(z) = G(z)z A z R (2.1.44) 

where the diagonal matrix A and a nilpotent matrix R are as above and the matrix valued 
function G(z) is analytic and invertible on a neighborhood of the point z — 



Exercise 2.1.14 Find a fundamental system of the form (2.1.44) for the system 



/ 



1 z 



zy = i o o ' y - 



Since the upper triangular matrix R is nilpotent the function 

^ = Id + i?logz + i? 2 ^ + ... 
is a polynomial in \ogz. We arrive at important 

Corollary 2.1.15 Solutions of a Fuchsian system have at most polynomial growth at z = 0. 



Definition 2.1.16 We say that an isolated singularity at z = of the system (2.1.1) is 
regular if all solutions to this system have at most polynomial growth at z — ► 0. 

From the previous results it follows that 

Fuchsian singularity =>■ regular singularity. 

The converse statement is false, as it follows from a simple 
Counterexample. The system 

*v=(? t) y (2 - L45) 
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has a singularity at z = of Poincare rank 1. However this is a regular singularity. Indeed, 
for the components y±, y 2 one obtains the following system 



y'l = 2/2 
z 2 y' 2 = Vi - zy 2 - 

The substitution y 2 = y[ into the second equation gives Euler equation for y\ , 

z 2 y'{ + zy[-yi = 0. 

Solving the latter 

yi = ci z H 

z 

yields 

C2 

z z 

So the fundamental matrix reads 

' z 1/z 



y(z) 



1 -l/z 



In general the following statement holds true. 



Theorem 2.1.17 A system of the form (2.1.1) with an isolated regular singularity at z = 
is weakly equivalent to a Fuchsian system. 

The proof of this theorem consists of two parts. First, we establish existence of a funda- 
mental matrix of the form 

Y(z) = G(z)z L 

with some matrix L and a matrix valued function G(z) analytic and invertible on a punctured 
disk (see Exercise 2.2.6 below). Since the singularity is regular we derive that G(z) has at 
most a pole at z = 0, 

G(z) = ±M(z) 

for some r > and a matrix M{z) analytic for \z\ < r and nondegenerate for \z\ > 0. For 
the sake of simplicity of notations let us do a shift 

L i-> L - r ■ Id 

in order to recast the fundamental matrix into the form 

Y(z) = M{z)z L . 

If det M(0) / Owe are done: the coefficient matrix 

A(z) = Y'(z)Y- 1 (z) = M^ M~ 1 + M' M- 1 

has a simple pole at z = 0. 

It remains to consider the case of degenerate matrix M(0) (recall that det M{z) ^ for 
\z\ > 0). The crucial point in the proof is the following statement, sometimes called Sauvage's 
Lemma. 
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Lemma 2.1.18 Given a matrix valued function M(z) analytic for \z\ < r, invertible for 
\z\ > such that det M(0) = 0. Then there exist a matrix P(z) polynomially depending on z 
satisfying 

detP(z) = l, (2.1.46) 

a matrix M(z) invertible and analytic on \z\ < r, and n nonnegative integers k±, . . . , k n such 
that 

M(z) = P(z)z K M(z), K = diag(k 1 ,...,k n ). (2.1.47) 

In order to complete the proof of the Theorem let us apply to the original system a gauge 
transformation 

Y(z) = P{z)z K Y(z). (2.1.48) 
The new system will have a fundamental matrix of the form 

Y(z) = M{z) z L with M(0) + 0. 

So the new system will have a Fuchsian singularity at z = 0. It remains to observe that 

det (P(z)z K ) =z fe i+-+ fe «. 

Therefore (2.1.48) is a weak equivalence of the original system with a Fuchsian system. □ 

The proof of this Lemma can be found, e.g., in §11 of Chapter IV of Hartman's book. 

Exercise 2.1.19 Find a weak equivalence of the system (2.1.45) with a system Fuchsian at 
z = 0. 

At the end of this section we consider the case of scalar linear differential equations of 
higher order 

+ ai(z)y( n -V + ■ ■ ■ + a n (z)y = (2.1.49) 
with coefficients analytic on a sufficiently small punctured disk 

< \z\ < p. 

Exercise 2.1.20 Let the coefficients satisfy the following conditions near z = 0: 

z z ai(z) is analytic at z = 0, z = l, 2,...,n. (2.1.50) 
Prove that the substitution 

yi = y 

V2 = z y'± 

(2.1.51) 

Vn = zy'n-l 
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reduces the equation (2.1.50) to a system with a Fuchsian singularity at z = of the form 



( o 




l 
o 





V -b n (z) -b n -i(z) 



... \ 

1 ... 

1 

... -b 2 (z) -h(z) ) 



with some functions b\(z), . . . , b n (z) analytic at z = 0. 



Hint: rewrite the equation (2.1.50) in the form 



d_ 

dz 



y + h(z) z 



A. 

dz 



y + ...b n (z)y = 0. 



(2.1.52) 



(2.1.53) 



Remark 2.1.21 The characteristic equation 

\ n + 6i(0)A n - 1 + • • • + b n (0) = (2.1.54) 

for the eigenvalues of the unperturbed matrix for (2.1.52) can be also obtained from the limiting 
Euler equation 

z n y^ + a^y^ + ■■■ + a n y = (2.1.55) 

where 

cti := lim z l ai(z), i = l,...,n. 



Definition 2.1.22 Under the assumptions (2.1.50) we say that the n-th order linear differ- 
ential equation has a Fuchsian singularity at z = 0. 

Differently from the case of systems the notions of regular and Fuchsian singularities 
coincide for scalar differential equations, as it follows from the following theorem, due to 
L.Fuchs. 

Theorem 2.1.23 All solutions to the scalar differential equation (2.1.49) with coefficients 
a±(z), a n (z) analytic on some punctured disk < \z\ < r have regular singularity at 
z = iff the singularity is Fuchsian. 

Proof: Sufficiency of the Fuchsian condition follows from the reduction of a Fuchsian scalar 
differential equation to a Fuchsian system given in Exercise 2.1.20, due to Corollary 2.1.15. 
Let us prove necessity. 

We begin from a simple lemma about multivalued functions analytic on the punctured 
disk. 



Lemma 2.1.24 Any linear differential equation (2.1.49) with coefficients analytic on the 
punctured disk < \z\ < r has a solution of the form yo(z) = g(z)z x for some complex A with 
a function g(z) analytic on the punctured disk. 
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See Lemma 2.2.5 below for the proof. 

As all solutions to the equation have regular singularity at z = 0, the function g(z) has 
at most a pole at this point, 

g{z) = -^jr for some k G Z, fr(z) is analytic for |z| < r, /i(0) / 0. 
Doing a shift A i— > A — k we obtain the solution in the form 

yo(z) = h(z)z x . 

We will now prove the Theorem using induction in n. For n = 1 the statement is obvious: 
the coefficient a{z) of a first order linear differential equation 

y + a{z)y = 

reads 

y' A fc'(z) 

a(z) = = — -. 

V z h(z) 

The second term in this formula is analytic at z = 0. 
To do the inductive step let us do a substitution 

V 

u = — 

yo 

to the equation (2.1.49). We obtain a linear differential equation for the function u that has 
a trivial solution u(z) = 1. Thus the coefficient of u vanishes in this new equation: 

u {n) + bi(z)u {n - 1 ^ + ■■■ + b n -i{z)u' = 0. (2.1.56) 

We arrive at a differential equation of order n — 1 for the function v = v! . Solutions to 
this equation have regular singularities at z = 0. From the inductive assumptions it follows 
existence of limits 

lim z % bi(z), i = 1, . . . , n — 1. 

Let us compare the coefficients bi{z) with the coefficients of the original system. To this 
end we plug y = uyo into (2.1.49) and collect the coefficients of u^ n \ u( n ~ l \ . . . , v! . After 
division by yo we obtain a triangular transformation 



i 



yo 



, i n \ y'o f n-l\ y' 



n-l J y \n-2 J y 
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where we have used the Leibnitz rule for derivating products 
(„„)(*) =u( fc W( \ ) n^V+j^ ^ ) u (fc - 2 V' + ---+( J ) uV*" 1 ) (2.1.57) 

Derivating the function yo = h(z)z :K , h(0) / it is easy to see that 

(k) 

yo 

has a pole of order k at z = 0. Since the functions &i, 62, . . . , & n -i have poles of the orders 
1, 2, . . . , n — 1 due to the inductive hypothesis, we easily conclude the same order of poles 
for the functions a±, 02, . . . , a n -\ respectively. Finally, from the equation for yo 

(n) (n-1) / 

h a\ !-■■• + a n -i \- a n — U 

yo yo yo 

we prove that a n has a pole of order at most n. □ 

Example. Bessel equation 

x 2 y" + xy +{x 2 - v 2 )y = 0, v G C. (2.1.58) 
The solution has a Fuchsian singularity at x = 0. Let us look for a solution in the form 

y = x A (c + cix + c 2 x 2 + ...). 

We have 

j/ = x x ~ 1 [Ac + (A + l)cix + (A + 2)c 2 x 2 + ...] . 

In a similar way 

y" = x A ~ 2 J^(A + fc)(A + fe — 1) c fc x fc . 

fc>0 

After the substitution and division by x A we obtain 



= 0. 



^ [(A + fc)(A + k - 1) c fc x fc + (A + fc) c k x k + (-i/ 2 c fc + c fc _ 2 ) 

fc>0 

At k = we have 

A(A-l) + A-z; 2 = 0, 

so 

A = ±v 2 . 

Assume that v 7^ (thus no Jordan blocks occur!). Take A = v. We obtain 

[{u + 1) 2 - v 2 ] Cl = o^ci = 0, 

[(1/ + A;) 2 - v 2 ] c fc + c fe -2 = 0. (2.1.59) 

So 

Codd = 
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if 2u is not a negative integer. Choose 

co 



2T(i/+ 1)" 
Then 

C2k ~ 



2 2k+v k\T{v + fc + 1) 
where we have used the recursion relation for the gamma-function 

r(i + i) = ir(s). (2.1.60) 

We obtain finally a solution to the Bessel equation in the form 

fc=0 v ; 

The series converges for all x. The sum of this series is called the Bessel function (of the first 
kind). If v Z then the functions J u {x) and J_^(x) give a fundamental system of solutions 
to Bessel equation. For u = one can obtain a second solution to Bessel equation linear 
independent with Jo(x) taking 

K (x) := ^J v (x)\ v=0 . (2.1.62) 

Exercise 2.1.25 Prove the following recursions for Bessel functions 

^-[x u J v (x)] = x v J v - 1 (x) (2.1.63) 
ax 

^ [x-"J v (x)] = -x' v J u+1 {x). (2.1.64) 



2.2 Monodromy of solutions to linear differential equations, local theory 

Recall (see Section 1.2 above) that monodromy describes permutations of branches of mul- 
tivalued analytic functions. In this Section we will study monodromy of functions defined 
on a punctured disk < \z\ < r for a sufficiently small r (that is, analytic functions on the 
universal covering of this punctured disk) . We begin with the following simple 

Example. Consider analytic function z x , A G C. From the definition 

z \ ._ e A(log|z|+iargz+27r in) n £ % 

it follows that branches of this function are labeled by integers. In order to determine the 
result of analytic continuation of this function let us consider a closed loop 

z(t) = ze 2wi \ 0<t<l 
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oriented counter-clockwise. We have 

^(•^jA _ e A(log|z|+iarg2:+27r it+2w in) ^ 2 1) 

At t = 1 we obtain the result of analytic continuation 

= gA(log|z|+iarg2:+27ri(n+l)) _ ^(0)] A^tt i A^ 

More generally we will adopt the following notations: the result of analytic continuation 
of a multivalued function f(z) analytic on the punctured disk we will denote 

For example, for f(z) = sj~z 

f(ze 2 ^)=-f(z); 

for f(z) = log z 

f{ze^ i )=f(z)+27ri. 

Let us study the monodromy of solutions to a system of n linear differential equations 

V = A(z)y (2.2.2) 

with the matrix of coefficients A{z) analytic on the punctured disk < \z\ < r. The following 
simple statement is the starting point for subsequent considerations. 

Lemma 2.2.1 Given a solution y(z) to the system (2.2.2), then the result y (ze 2ir! ) of an- 
alytic continuation of the solution along the loop (2.2.1) is again a solution to the same 
system. 

Proof: After the substitution z = ze 2 " we obtain for the function y := y(z) the system 

Because of analyticity of A we have A(z) = A(z). □ 

Definition 2.2.2 The operator of analytic continuation 

y(z)^y{ze 2 * 1 ) (2.2.3) 

acting on the space of solutions to the system of linear differential equations (2.2.2) is called 
the monodromy operator of (2.2.2) 

As the space of solutions to the system of n linear differential equations is n-dimensional, 
the monodromy operator is represented by a n x n matrix M in any basis of solutions. Putting 
the basic solutions in the columns of a fundamental matrix Y{z) one spells out the definition 
of the monodromy matrix in the form 

Y (ze 2wi ) =Y{z)M (2.2.4) 
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Because of nondegeneracy of the fundamental matrix the monodromy matrix does not de- 
generate. 

Example. The monodromy of the fundamental matrix (2.1.36) around a Fuchsian sin- 
gularity is a diagonal matrix 

Y(z) = G{z)z A ^Y (ze 2lTi ) = Y{z)M, M = e 2ni A = diag (e 2ni A \ . . . , e 2wi A ") . (2.2.5) 



Exercise 2.2.3 Prove that the monodromy of the fundamental matrix (2.1.44) is described 
by the matrix 

M = e 2wiA e 2wiR . (2.2.6) 
Hint: Check that the matrices e 2nlA and e 2nlR commute. 

Observe that the matrix e 2ntR is upper triangular; all its diagonal entries are equal to 1. 

Exercise 2.2.4 Given a non- degenerate complex nxn matrix M, prove existence of a matrix 
L such that 

M = e 27riL . (2.2.7) 



The following statement is often useful in the study of local properties of solutions to 
linear differential equations near a singular point. 



Lemma 2.2.5 There always exists a solution to the linear system (2.2.2) near an isolated 
singularity z = that can be represented in the form 

y(z)=g(z)z x (2.2.8) 

for some A G C where the vector valued function g(z) is analytic on the punctured disk 
< \z\ < r. If the solution y(z) grows at most polynomially at z — ► then the vector valued 
function g(z) has at most a pole at z = 0. 



Proof: Let y{z) be an eigenvector of the monodromy operator with an eigenvalue ji. Because 
of nondegeneracy of the monodromy matrix ^ / 0. Put 

A = ^— log n 

for some choice of a branch of the logarithm. The product 

g(z) := y{z)z~ x 

is a single-valued analytic function on the punctured disk. It remains to recall that the mon- 
odromy operator on the finite dimensional space of solutions always possess an eigenvector. 

In order to prove the second part of Lemma it suffices to observe that the function g{z) 
grows at most polynomially at z — > if the function y(z) does so. Since g(z) is analytic on 
the punctured disk it must have a pole at z = 0. □ 
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Exercise 2.2.6 Prove that a system of the form (2.2.2) with an isolated singularity at the 
origin possesses a fundamental matrix of the form 



Y(z) = G(z) z 1 



(2.2.9) 



for some constant matrix L and a matrix valued function G(z) analytic and nondegenerate 
on the punctured disk < \z\ < r. Prove that in the case of regular singularity the entries of 
the matrix G(z) have at most poles at the origin. Deduce that coefficient matrix A(z) must 
have at most a pole at the origin if the singularity is regular. 



2.3 Monodromy of Fuchsian systems and Fuchsian differential equations 

We now proceed to the global study of monodromy of differential equations with rational 
coefficients. 

A system of n first order linear differential equations with rational coefficients 

y' = A(z)y (2.3.1) 

is said to be Fuchsian if it can be represented in the form 

(z-z )y'= [A + A t {z - z ) + A 2 (z - z ) 2 + . . .] V (2.3.2) 

near any singular point z$ £ C. 

Equivalently, using the simple fraction decomposition of the rational matrix valued func- 
tion A(z) one obtains a representation of a Fuchsian system in the form 

dy 



dz 



Ai A k 
_J_ + ... + *_ 



z — Z\ 



z- z k 



y 



(2.3.3) 



where Ai , . . . , A n are n x n constant matrices and the complex numbers z± , . . . , z k are 
pairwise distinct. The singularities of the system (2.3.3) are at these points and, possibly, at 
infinity. 



Lemma 2.3.1 After the substitution 
the system (2.3.3) rewrites 



w = 



dy 
dw 



A r 



w 



0(1) 



where 



Aoo := -(Ai + ■■■ + A k ). 



(2.3.4) 
(2.3.5) 



Proof: After the substitution one obtains 

2 dy Ai 
—w — = w 
dw 



w 



1 — Z\W 

A\ H h A k + w 



+ ■■■ + 



A k 



1 - z k w 
z\A x 



1 — Z\W 



+ ••• + 



ZkA k 
1 - z k w 



y- 



(2.3.6) 
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So 



dy 
dw 



'A r . 



w 



ziAi 

1 — Z\W 



+ ... + 



z k A k 
1 - z k w 



□ 



Corollary 2.3.2 The infinite point of the Riemann sphere is a Fuchsian singularity of (2.3.3) 
if / 0; the system (2.3.3) has no singularity at infinity iff = 0. 



Exercise 2.3.3 Prove that after an arbitrary fractional linear transformation 

az + b 

w = :, ad — bcy^0 

cz + a 

the Fuchsian system (2.3.3) goes to another Fuchsian system with the singularities at 

az i + h ■ n , , n i 

Wi = :, t = 1, . . . , k + 1, where z k +i = oo 

czi + a 

and with the residue matrices A±, . . . , A^, A^+i := A^. 

We are ready to define an important characteristics of a Fuchsian system: the monodromy 
representation. Let us fix a nonsingular point 

zo G C \ {zi, . . .,z k ,z k +i}. 

For any closed loop 

7 : [0, 1] -> C \ {zi, ...,z k , z k+1 }, 7(0) = 7(1) = z 



denote M 7 the linear transformation 



V » M 7 (y) (2.3.7) 



the result of analytic continuation of a solution to (2.3.3) along the loop 7. According to 
Lemma 2.2.1 the result of the analytic continuation defines a linear operator in the space of 
solutions to (2.3.3) that does not depend on the deformation of the closed loop keeping the 
end points 7(0) = 7(1) = zq fixed. We obtain a map 

of the fundamental group 

TTi (C \ {zi, . . . , z k , z k+1 }; zq) -► Aut(C n ) (2.3.8) 

to the group of automorphisms of the linear space of solutions to the system (2.3.3). 

Choosing a basis in the space of solutions one can describe the monodromy (2.3.8) by 
matrices. Namely, analytic continuation of a fundamental matrix Y(z) along the loop 7 gives 



62 



a new fundamental matrix M 7 (Y) of the same system. The latter is related to Y by the right 
multiplication by an invertible matrix M 7 : 

M 7 (y(z)) = Y{z)M 1 . (2.3.9) 

Recall that the elements of the fundamental group tt\{X]Zq) of a topological space X 
with respect to the marked point zq £ X are the homotopy classes of the continuous loops 

7 :[0,1]-X, 7(0) =7(1)- 

The product 7 = 7172 of two elements in the fundamental group is defined by the consecutive 
running along the two loops 71 (i) and 72 (i): 

' 71 (2 *), < t < i 
7 (t) = { (2.3.10) 
, 72(2t-l), § <t < 1 

for < t < 1. The inverse element to 7 is defined as the homotopy class of the same loop 
passed in the opposite direction 

r\t) := 7 (l-t). (2-3.11) 
The neutral element is the trivial loop 

7 (i) = zo for any < t < 1. 

We are ready to prove the main property of the monodromy map (2.3.8). 

Lemma 2.3.4 The monodromy map (2.3.8) is an antihomomorphism of the fundamental 
group of the punctured Riemann sphere to the group of automorphisms of the space of solutions 
to (2.3.3), i.e., for any two loops 71, 72 on C \ {z±, . . . , Zk, Zk+i} one has 

M 7172 = M 72 M 71 . (2.3.12) 

Moreover, for any loop 7 

M 1 -x = M~ x . (2.3.13) 

To the homotopy class of the trivial loop it corresponds the identity map of the space of 
solutions. 

Proof: It suffices to prove validity of (2.3.12). Let us do it using the matrix realization 
(2.3.9) of the monodromy map with respect to a choice of the fundamental matrix Y{z). The 
analytic continuation of Y{z) along the loop 71 transforms Y(z) to Y{z)M~ fl . The subsequent 
analytic continuation of the latter matrix along 72 gives 

(Y(z)M l2 ) M 11 =Y{z)M l2 M 11 . 

This proves that 



M 7172 = M 72 M 71 . 



□ 
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Definition 2.3.5 Let Y(z) be a fundamental matrix of the system (2.3.3). The antihomo- 
morphism 

7ri (C\{zi,...,z k ,z k+1 };zo) -> GL(n;C) 

(2.3.14) 
(2.3.15) 

mapping the homotopy class of aloop 7 to the monodromy matrix M 7 is called the monodromy 
representation of the system associated with the chosen fundamental matrix. 

Exercise 2.3.6 Prove that a change of the fundamental matrix 

Y{z) 1 ^ Y(z) = Y{z)C 
changes the monodromy representation to the conjugate one: 

M 1 = C' l M 1 C. (2.3.16) 

Scalar case. Fuchsian differential equations of order n can be described according to the 
following statement. 

Lemma 2.3.7 Any Fuchsian differential equation of order n must have the form 
y^ + a 1 (z)y^ + --- + a n (z)y = 

(2.3.17) 

/ \ Qm\Z) 

a m(z) = m , m = l, 2,...,n 

P(z) = (z - zi)(z - z 2 ) ■ ■ ■ (z - z k ), Zi^Zi for i^j 

Qm(z) is a polynomial of degree m (k — 1), m = 1, 2, . . . , n. 
TTie regular singularities are at the points Z\, . . . , z k , z k+ \ = 00. 

Denote 

the characteristic exponents of (2.3.17) at the singular point z = Zi, i = 1, . . . , k, k + 1. 

Exercise 2.3.8 (c/. Exercise 2.3.3). Prove that after a fractional linear transformation 

~ az + b A u / n 
z^z = _, aa — bcf^u 

cz + a 

the Fuchsian equation (2.3.17) transforms to another Fuchsian equation with singularities at 

z% = — , , i = 1 , • • ■ , k + 1 
czi + a 

with the same characteristic exponents. 
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Theorem 2.3.9 (Fuchs relation). 
equation (2.3.17) satisfy 



The characteristic exponents of a Fuchsian differential 



k+1 

£ 

i=i 



n(n — 1) 



(fc-1) 



(2.3.18) 



Proof: For any system yi(z), . . . , y n (z) of n functions denote 



y[ 



W(z) = det 



y'l 



\y\ 



(n-l) 



V2 

y'2 
y'i 



(n-l) 
^2 



(n-l) 



(2.3.19) 



the Wronskian determinant of this system. The following properties of the Wronskian can 
be easily established: 

• If yi(z), . . . , y n {z) is a basis in the space of solutions to an n-th order linear differential 
equation (2.3.17) then the Wronskian W{z) never vanishes for z ^ z±, . . . , Zk, Zfc+i- 

• Given two bases yi(z), . . . , y n (z) and yi(z), . . . , y n (z), then the two Wronskians W(z) 
and VF(z) are related by multiplication by a nonzero constant, 

W(z) = pW(z), p^O. 



• The analytic continuation of solutions yi(z), . . . , y n (z) along a closed loop 7 G C \ 
{z±, . . . , Zk, ^fc+i}) 7(0) = 7(1) = zq results in multiplication of the Wronskian by a nonzero 
constant 

W(z ) - ^W(z ) 
independent on the choice of the basepoint zq. 

z = f(z) yields the following transformation 



• Change of the independent variable z 
of the Wronskian: 



W{z) 



n(n — l) 
2 



W(z). 



(2.3.20) 



Indeed, one has 



dy dz dy 
dz dz dz 



fy 

dz 2 



d 2 y d 2 z dy 
dz 2 dz 2 dz 



etc., and, more generally, 
d m y 



dz r ' 



dzy 1 d m y 
dz 



dy 



, h linear combination of 

dz m dz 



d 



m—l„ 



dz m - 1 ' 



m 



1, ... ,n - 1. 



So, up to adding of a linear combination of the previous rows, the rows of the Wronskian will 



be multiplied by 1, , (f|) 2 • • • , (gf )™ 1 resp. This gives the transformation law (2.3.20). 



i;\n-l 
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We will now study the behavior of the Wronskian at the singularities. Let us assume for 
simplicity that all characteristic exponents are pairwise distinct and nonresonant. Denote 

yf = gf{z){z - zi)^\. . .,y® = gM(z)(z - Zl ) X " , i = 1, . . . , k, k + 1 (2.3.21) 

the bases in the space of solutions associated with these exponents, i.e., the functions gf\z), 
. . . , gn\z) are analytic near z = zi and, moreover, gm(zi) = 1 for any m = 1, . . . , n. Without 
loss of generality we may assume, doing if necessary a fractional linear transformation, that 
all the k + 1 singular points are away from infinity. Then it is easy to see that the Wronskian, 
calculated for the basis y^\z), . . . , yn\z) behaves like 

W(z) =c l (z-z l ) x ^ + - +x -- 2i ^ 11 (l + 0(z- Zi )), z^ Zl . (2.3.22) 

where 

« = n ( a ? - a ?) + o- 

p<q 

At infinity, according to (2.3.20), one has 

W(z) = Coo z n ^ (l + O^yj, Coo + (2.3.23) 

(do a transformation 5 = 1/z). 

We now consider the logarithmic differential 

uj := dlogW(z). (2.3.24) 

The differential u is meromorphic on the Riemann sphere. It has poles at the points z\, . . . , 
Zk+i with residues 

veSz=ZiU = x f + ... + \f- n ^^, « = 1 fe + 1 

and infinity, 

res 2=00 u; = — n(n — 1). 

Applying the residue theorem 

resp u; = 

Pec 

one obtains (2.3.18). □ 



Exercise 2.3.10 Assume as above zt+i = oo. Prove that the substitution 

y = (z-z 1 ) ai (z-z 2 ) a2 ...(z-z k ) a *y 

transforms the Fuchsian equation (2.3.17) into another Fuchsian equation for y = y(z) with 
the same singularities and with the characteristic exponents 

= A m " «i> m=l,...,n, i = l,...,k. 
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Let us consider in details the case of second order Fuchsian equations with three singu- 
larities at z = zi, z = Z2, z = Z3. The solutions y = y(z) to this equation are conveniently 
denoted using the so-called Riemann scheme 





( 


Zl 


Z2 


Z3 


\ 


y = V 






a! 2) 


\(3) 
A l 


z 




\ 


4" 


X (2) 
A 2 


X (3) 
A 2 





where we keep the same notations as above for the characteristic exponents. To motivate 
this notation we will show that the equation is uniquely determined by the characteristic 
exponents. So the symbol (2.3.25) denotes a two-dimensional space of solutions to the asso- 
ciated differential equation. Recall that, due to Fuchs relation the characteristic exponents 
are constrained by the equation 

A « + A « + A W + A W + A (3) + A (3) = L {2 3 26) 



Theorem 2.3.11 (Riemann). The coefficients of the second order Fuchsian equation with 
three regular singularities are uniquely determined by the position of singularities and by the 
characteristic exponents. 

Proof: The equation must have the form 

y" + p(z)y' + q(z)y = 0. (2.3.27) 

ai a 2 a 3 
p{z) = + + 



q(z) 



Z — Zl Z — Z2 Z — Z3 

1 / bi , 62 &3 



(z — Zl)(z — Z2)Z — Z3) \Z — Zl Z — Z2 Z — Z3 



a J = l-(Af ) + A«) 



bi = X^X 2 l \zi - Zj)(zi - z k ), i = 1,2, 3. 



□ 



Exercise 2.3.12 Prove that a second order Fuchsian differential equation with singularities 

(i) (i) 

at zi, . . . , Zk, Zk+i = 00 with the prescribed characteristic exponents X\ , X 2 , i = 1, . . . , k + 1 



67 



must have the form 

y" + p(z) y' + q(z) y = 



(2.3.28) 



= E — TTT^ 



g(z) = 



i=i 
1 



V p^.) + A (o°) A (oo) ^-2 + + . . . + ^ 



i=l 



P(z) 

where P(z) = JJ(z — Zj). 
i=i 

Generalize to the case of Fuchsian equations of higher orders. 



For > 2 the coefficients ci, . . . , c^_2 are not determined by the poles and the exponents. 
They are called accessory parameters of the Fuchsian equation. 



2.4 Gauss equation and hypergeometric function 

x(l - x) y" + [7 - (a + + 1) x] y -af3\ 



= 0. 



Riemann scheme for the Gauss equation: 

y = V 



1 00 

a ; x 

1 — 7 7 — a — /3 /3 



(2.4.1) 



(2.4.2) 



Let 
Denote 

i.e. 



7/0, -1, -2,... 
T(a + n) 



(a) n := 



r(«) ' 



(2.4.3) 



(a) = 1, (a) n = a(a + 1) . . . (a + n - 1), n = 1, 2, 
Define the hypergeometric series by 

" («)n(/?)n x n 



2 F 1 (a,P;r,x) := ^2 



n=0 



(2.4.4) 



Lemma 2.4.1 The series (2.4.4) converges for \x\ < 1. /t gives i/ie unique solution to Gauss 
equation (2.4.1) analytic at x = 0. 
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The resulting function 2-^1(0, P\ 7! x) analytic for |x| < 1 here will be also denoted 
F(a, (3; 7; x) for brevity. 

The hypergeometric function can be analytically continued to the entire complex plane 
with a branch cut from x = 1 to x = 00. This gives a possibility to construct a basis of 
solutions to the hypergeometric equation near x = 1 and x = 00. 

The substitution 

y = x 1_7 -u 

yields a Gauss equation for u of the form 

x(l - x) u" + [2 - 7 + (a + /? - 2 7 + 3) x] v! - (a - 7 + l)((3 - 7 + 1) u = 0. (2.4.5) 
So the function 

y = a; 1 -')' J p(a-7 + l,/3-7 + l;2-7; x), 7 0Z (2.4.6) 
gives another solution to (2.4.1) corresponding to the characteristic exponent 1 — 7. 
In a similar way at x = 1 one obtains the following two solutions 

y = C 1 F(a,0;a+p—y+l; 1 -x) +C 2 (1 -xf-^F^-a, 7-/?; 7 -a-/?+l; 1-x). (2.4.7) 
At infinity: 

y = Cix~ Q F(a, 1 - 7 + a; 1 - /? + a; x" 1 ) + C 2 x^F(p, 1 - 7 + /3; 1 - a + /?; x" 1 ) (2.4.8) 

Exercise 2.4.2 Derive the following particular cases of the hypergeometric function: 

F(-n, P; p; -x) = (l + x) n , n£Z (2.4.9) 

F(l,l;2;-X) = l0g(1 x +X) . (2.4.10) 



Exercise 2.4.3 Prove that 



-^F(a,P; T ,x) = — F(a + l,/3+l; 7 + l;x). (2.4.11) 

ax 7 



Exercise 2.4.4 For the case 7 = 1 prove that the second solution of Gauss equation linearly 
independent from F(a, P; 1; x) reads 

y = lim( 7 - l)- 1 [x x -TF(a - 7 + 1, P - 7 + 1; 2 - 7; x) - F(a, p; 7; x)] 
7— >i 

(2.4.12) 

^, ftl:I >>o g , + £<^»£(^ + ^-?). 

n=l fc=l x 7 
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Monodromy of Gauss equation, "brute force" method: look for three matrices Mq, Mi, 
Moo with the eigenvalues (l,v~ l ), (l,z/A~V *) an d (A, fi) respectively satisfying 



M OD M l M = Id. 



Here 



ia 



A = e 

The nonresonancy will be assumed: 

A^/i, v 7^ 1, v 7^ Xfi. 
Without loss of generality we may assume M^ to be diagonal, 



Denoting 



Mi 



a b 
c do 



A 
n 



Mi 



ai 6i 



one obtains a system of equations 



A (a ai + &1C0) = 1 
^i(6 ci + dorfi) = 1 
ai&o + hdo = 
aoci + c d\ = 0. 



a 



1 



a Q d - b c = v 1 



«i + &i = 1 + v A V 1 , 



aidi — 6iCi = ^A 1 /i 1 . 



After some calculations one obtains 

/ M (A _^_ 1) + 1 



Mi = 



A-/i (A-IKm-IKA-z^-i/) A 
\ ~ Aiz(A-Ai) 

(A-l)( M -l)(^- M )(A-z/) 



X — fi 



^(M-A) \ 
^( I/ _ Ai + l)_l y 

A-/x \ 



V 



A-I/+ £ 



(2.4.13) 
(2.4.14) 



(2.4.15) 



(2.4.16) 



Observe that the monodromy matrices Mo, M^ are determined non uniquely: there remains 
an ambiguity doing simultaneous conjugations by a diagonal matrix D: 

(M , Mi) i ^ (D^MoD, D~ 1 M\D) , D = diag(di, d 2 )- 

Such a conjugation does not change the diagonal matrix M^. 

Second method: using Barnes integral representation of hypergeometric functions. 
Consider the following integral 



r( 7 + s) 



T(-s)(-z) s ds. 



(2.4.17) 
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It is assumed that 

I arg(-z)| < vr; 

the integration path is chosen in such a way that the poles of the function T(a + s)r(/3 + s) 
are on the left from the integration path and poles of the function r(— s) are on the right 
from the integration path. Recall that the function 



r(ar) = / 
Jo 



oo 

e~H x ~ x dt, Rex>0 



can be analytically continued to a meromorphic function on the complex plane by using the 
functional equation 

r(x)r(i -x) = (2.4.18) 

sm TTX 

It has poles at the non positive integer points 

x = 0, -1, -2, ... 

with the residues 

ie Sx= . n T(x) = (2.4.19) 
For large |x| the asymptotic behaviour of gamma- function is described by the Stirling formula 

log T(x + a) = ^x + a—^j log x — x + ^ log 2ir + o(l) 

(2.4.20) 

for |x| — ► oo, | arg(x + a)\ < tt — e and | argx| < ir — e 

for any small positive e. 

From the Stirling formula it follows that, for large |s| on the integration contour the 
integrand behaves like 



O 



\s\ a+l3 7 1 exp{— arg(— z) • Ims — tt |Ims|} 



Therefore the integral converges and, moreover it defines a function analytic in z on the 
domain | arg z\ < ir — e. 

Let us prove that the function / = f(a, /3; 7; z) satisfies the Gauss equation in z. Using 
the identity 

xT(x) =T(x + l) 
we find for the derivatives of the function / the following expressions 

27riJ_ loo r( 7 + s ) v 
27TiJ_ ioo r( 7 + s ) k 
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After the substitution into the Gauss equation one obtains 
z(l - z) f" + [7 - (a + + 1) z] f - a/3f 



(2.4.21) 



2tt 

In the first integral do the substitution 

s = l + t 

The new integrand reads 

T(a + S )T((3 + s) (7 + a _ 1)r(1 _ fl)( _^-i 



r( 7 + s ) 

= r( ° + r ( ;^fi? +1) h + m ~ t){ - z) ' - ta + twi + " r(a ^f) + " r ""- t) ' 

So the two integrals in (2.4.21) cancel. 

We now want to prove that the Barnes integral is equal to the hypergeometric function 
(2.4.4), up to a constant factor. 

Lemma 2.4.5 For \z\ < 1 one has the following identity 

f(a, [3; 7; z) = r( ^ (/3) F(a, (3; 7; z) . (2.4.22) 



Proof: Let us consider the integral of the same expression as in (2.4.17) taken along the half 
circle C of radius N + \ , N € Z with the center at the origin laying to the right of imaginary 
axes. Using the functional equation (2.4.18) we arrive at the following integral 

J_ f ''•'V/^" (-,)•*. (24.23) 

From the Stirling formula we obtain that 

7TT(a + s)T((3 + s) . ( a+ Q_~i\ {—z) s 



r(7 + s)T(l + s) sin ns V / sin7rs 

uniformly in arg s on the contour C when N — > 00. The substitution 
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yields, for \z\ < 1, 

(-zy 



o 



siri7rs 



1 



exp<^ [ N + - cos8 log\z\ - [ N + - ) sin6> aig(-z) - [N + - tt | sin 0\ 



1 



= O 



1 



exp <M iV + - cos 6» log |z| - iV + - e | sin 6»| 



1 



O 



O 



exp{2-i {N + \) log|z|}] , < |^| < ^7T 
exp{-2-ie (AT + i)}] , Jtt < |0| < |tt. 



So, for log |z| < (i.e., for |z| < 1) the integrand goes to zero rapidly enough on the arc C. 
Therefore 

/ for iV^oo. 

Jc 

Let us now consider the following loop integral: 



/«oo /■— + 5 1 f moo 

" / +/ +/ 

-too J-ioo ic Ji(N+h 



("+*)J 

By Cauchy theorem it is equal to the sum of residues of the integrand at the poles s = 
0, 1, 2, . . . , N. Let iV go to infinity. Then the last three integrals will go to zero under 
assumption 

| arg(— z)\ < ir — e and \z\ < 1. 
It remains to compute the residue 

for any nonnegative integer n. Applying (2.4.19) we obtain 

T(a + s)T((3 + s) r/ T{a + n)T{f3 + n) z n 

res s=ra '-f '- r (-a) = — v ^ — . 

I (7 + s) T(7 + n) n! 

Finally we obtain for the Barnes integral the hypergeometric series 

, . A T(a + n)T(p + n) z n Via) Tiff) n 

Af^no^^ T(7 + n) n\ T(7) 



ra=0 



□ 



We are now ready to derive the formula expressing the hypergeometric solution F(a, (3; 7; z) 
with two solutions (2.4.8) defined near z = 00. Consider the integral 



1 /" T{a + s)T{f3 + s) 
2ni] D r( 7 + s) 



r(-s)(-z) s ds (2.4.24) 



where D is the half of a circle of the radius R with the centre at the origin lying on the 
left from the imaginary axis. Like in the previous calculation one can show that the integral 
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(2.4.24) goes to zero when R — > oo and \z\ > 1, assuming that | arg(— z)\ < it and the 
sequence of radii R is chosen in such a way that the distance from the arc D to the poles of 
the integrand is bounded from below by a positive constant. Computing the residues of the 
integrand at the poles 

s = —a — n and s = — (3 — n, n G Z 
and applying Cauchy theorem one finds 



i_ i*~ r ( q + .)r(g + .) 



oo 



Er(a + n)T(l — 7 + a + n) sin7r(7 — a — n; _ n _ 
r(l + n)T( 1 - + a + n) cos vrn sin ir(3 - a - n) ' ^ 



n=0 



"(1 + n)T(l — /3 + a + n) cos 7rn sin ir((3 — a — n) 



rQg + n)T(l - 7 + (3 + n) sin7r( 7 - /? - n) , w/3-n 
^ T(l + n)r(l - a + /3 + n) cos7rn sin7r(a - (3 - n) Z 



n=0 



Finally one arrives at the following connection formula expressing one solution to Gauss 
equation as alinear combination of two other solutions: 

y F(a,(3;r,z) = J, (-*) -F a, 1 - 7 + a; 1 - /3 + a; 



r( 7 ) r(a- 7 ) v 7 V ' ' ' ' 'z, 

(2.4.25) 

The expression (2.4.25) is valid under the assumption 

I arg(— z)\ < 7r. 

Exercise 2.4.6 Compute the integral 

1 pi 00 

-j_j (a + a)n - a){ - zfda 

assuming | arg(— z)\ < n. 

Exercise 2.4.7 (Gauss formula). Prove that 

/or 

Re (7 - a - (3) > 
Hint: use Euler integral representation (2.4.37). 
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Exercise 2.4.8 (Barnes Lemma). Prove that 

I pi oo 

/ r( a + s)T((3 + s)T(j-s)T(S-s)ds 

27TlJ_ ioo 

(2.4.27) 

_ F(a + 7) Tja + 5) T{P + 7) T(J3 + 5) 
T(a + P + j + 5) 

The integration contour is chosen in such a way that poles o/T(7 — s) T(5 — s) are on the left 
and poles ofT{a + s)T{(3 + s) are on the right of it. (It is also assumed that the nmbers a, 
(3, 7, 5 are such that neither of the poles of the first function coincides with any pole of the 
second function. 

Exercise 2.4.9 Using Barnes Lemma prove validity of the following connection formula 
T( 7 - a) T( 7 - (3) T(a) T(0) F(a, (3; 7; z) = 

(2.4.28) 

= T( 1 )T(a)T((3)T( 1 -a-P)F(a,P;a + P- 1 + l;l-z) + 

+r( 7 ) T( 7 - a) T( 7 -P)T(a + P- 7) (1 - z) 7 " a - /3 F( 7 - a, 7 - /?; 7 - a - (3 + 1; 1 - z) 
assuming that 

I arg(l — z)\ < 7T and |1 — z\ < 1. 

Hint: 

We finally arrive at the following result. Assume that none of the numbers 1 — 7, 7 — a — /3, 
a — j3 is an integer. Then we can construct three bases in the space of solutions to Gauss 
equation as it was explained above (see formulae (2.4.6) - (2.4.8)). We will write the three 
pairs of basic vectors as three row matrices 

Y (x) := (F(a, (3; 7; x), x^Fipt - 7 + 1, /3 - 7 + 1; 2 - 7; x)) (2.4.29) 

Yl(x) := (F(a,/3;a + /?-7 + l;l -x),(l - x^^F^ - a, 7 - /3; 7 - a - /3 + 1; 1 - x)) 

(2.4.30) 

Yoo(x) = (x- a F{a, 1 - 7 + a; 1 - + a; x" 1 ), x^F(/3, 1 - 7 + 0; 1 - a + /?; x" 1 )) (2.4.31) 

In these formula we assume that x belongs to the upper half plane Imx > 0. Then in the 
definition of fractional powers of x and 1 — x we choose the principal branch of the arguments 

< argx < 7r, < arg(l — x) < it. 
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Theorem 2.4.10 Under the above assumptions the following connection formulae hold true 
for Im x > 

r( 7 )r(7- a -/3) r(2- 7 )r(7-g-/3) \ 

r( 7 -a)r( 7 -/3) r(l-a)r(l-/3) 

y (x) =yi(x) | (2.4.32) 

r( 7 ) r(a+/3- 7 ) r(2- 7 )r(a+/3-7) , 

r(a)r(/3) r(a- 7 +l) r(/3-7+l) / 



-i7ra r(7)r(^-Q) 
r( 7 -a)r(/3) 



i7r ( a - 7+ l) r(2- 7 )r(/3-a) \ 
e r(/3- 7 +l)r(l-a) \ 



-i7r/3 r( 7 )r(a-/3) i7r(/? - 7+ i) r(2- 7 )r(a-/3) 



r( 7 -/3)r(a) 



r(a- 7 +i)r(i-/3) / 



(2.4.33) 



The formula (2.4.32) remains valid for |arg(— x)| < 7r; the formula (2.4.33) is valid for 
| arg(l — x)| < 7r. 

These connection formulae yield the following expressions for the monodromy matrices in 
the basis Y^: 

e -2iwj 
M Q = ^—^ X 



g2i7TQ g2i7r/3 



/ p 2iir(p+a) _ „2i7r(/3+ 7 ) _ p 2inf3 , _2«r7 47rV^ +^r(/3-q) \ 

/ e e e + e r(/3)r(i-a)r(a -/3)r(/3- 7 +i)r( 7 -a) 

47rV 7r ( 2a +T>r(a-/3) „2ivr(a+ 7 ) _ „2i7r(a+/3) i „2ina _ e 2in~/ . 

\ r( a )r(i-/3)r(^-«)r(a-7+i)r( 7 -/3) e e + e e / 



(2.4.34) 



M 



g2i7ra g2i7r/3 

^ g-2i7ra ^g2i7r(a+ 7 ) _ g2i7r(a+/3) _|_ R 2i-na _ e 2iir"/^ 



An e' 7r " l T(/3— a) 



V 



47r 2 e' 7r Tr(a-/3) 

T(a)r(l-/3)r(/3-a)r(a- 7 +l) r( 7 -/9) 



(2.4.35) 



r(/3)r(l-a)r(a-/3)r(/3- 7 +l)r(7-Q) 
e -2i7r/3 ^ e 2i7r(/3+a) _ e 2i7r(/3+7) _ e ^/3 + & 2i-K^ J 



„2ina 



Moo = 



e 2 ^ 



(2.4.36) 



The third method of computation of the monodromy of Gauss equation is based on the 
Euler integral representation of the hypergeometric functions: 



F(a,P;r,x) 



r(7) 



r(/3)r( 7 



— - / ^-^l-tJT-^-^l-tx) -01 ^, Re7>Re/3>0, |x| < 1. 

(2.4.37) 



Exercise 2.4.11 Prove (2.4.37). 
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Hint: expand (1 — tx) a in geometric series and use Euler formula 

T(x + y 

for beta-integrals 

B(x,y):= f t x ' l (l - t) y - l dt, Rea;>0, Rey > 0. 
Jo 

We will consider here only a very particular case of this representation. 

Exercise 2.4.12 Consider the complete elliptic integrals 

K= f 1 - dS (2.4.38) 

Jo s 2 ){l-k 2 s 2 ) 

-i/k 

^(l-s 2 )(l-fc 2 s 2 ) 
as functions of the complex variable 

x = k 2 , x 0, 1, oo. 

Prove that the functions y\ = K(x), y2 = K'(x) are two linearly independent solutions to the 
following Gauss equation 

x(l - x) y" + (1 - 2 x) y' - - y = 0. (2.4.40) 
Prove that in the basis y\, y2 the monodromy matrices have the following form: 



/■ i//<: (is 
iK'= = (2.4.39) 

Ji J(l- s 2 )(l-k 2 s 2 ) 



Derive the formula 



K = -F(l/2,l/2;l;k 2 ). (2.4.42) 
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The following books were useful for preparing the present lecture notes. 
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